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*•  ABSTRACT 

The    sputtering   process    has    been   investigated   by    simulating 

the    sputtering    of    single-crystal    copper    with    1-7    keV    argon.       A 

digital    computer    was    used    to   build    the    crystal,    bombard   it,    and 

** '    •  M I  --#W 

move    crystal   atoms.       Four    mechanisms    were    observed   which 
cause    surface   atoms    to    sputter.      An   atom    is    sputtered    when 
(1)    it   is    squeezed    out   of    the    surface,     (2)    it   is    scooped    out 
when   another   atom    strikes    its    inner   hemisphere,     (3)    it   is    ejected 
when   an   atom    passes   behind   it,    and    (1+)    it   is   knocked   out   by   a 
second   layer    atom    which    is    moving   outward.       Nearly   all    sputtered 
atoms    were    surface   atoms.       Second   and    third   layer   atoms    were 
sputtered   only   for    ion   energies    greater    than    5   keV .       They   were 
sputtered  by    mechanisms    similar    to    the    surface   atom    mechanisms. 
"Silsbee    chains"    were    observed    to   be    directed    into    the    crystal, 
and    momentum   focusing  was    observed   to    cause    sputtering   only 
when   it   occurred   in    close   packed,    surface   rows.      Outward 

directed    chains    were   not    observed.       Sputtering   deposit   patterns, 

i 

sputtering   ratios,    and    sputtered   atom    energy   distributions    were 

obtained   for    (100),     (110),    and    (111)    surfaces.      All   data 
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1.       INTRODUCTION. 

Material   is    lost   from    a    substance   undergoing   ionic   bombard- 
ment.      This    is    sputtering.      Many    authors   have    tried    to    define 
sputtering   using   terms    such   as    "( atoms )...  ejected   or   knocked 
out",    "emission",    "ionic    erosion",    "disintegration",    and   "break- 
down."     A   reader    who   blithely   accepts    one   of    these   definitions 
in   the    course   of   his    reading    may   conclude:       This   definition 
describes    the    mechanism    of    sputtering.       This    is    an   illogical,    if 
not   erroneous,    conclusion;    the    definition   of    sputtering   must  be 
derived   from    the    mechanism,    not   the    mechanism   from    the 
definition.       The    sputtering   process   has   been  investigated 
utilizing   one   of    two    general    mechanisms,    ejection   or    emission. 
Ejection   implies    an   immediate    or    direct   reaction   release   of    an 
atom    from    a    crystal;    emission,    a   delayed    or    indirect   reaction 
process    which   results   in   the   release   of    an   atom. 

In    1923,    Kingdon   and   Langmuir  bombarded   thoriated 

tungsten    with   various    ions    in   a    glow    discharge   tube.       This    was 
a    special    case   of    sputtering   since    the    thin   surface   film    of 
thorium    on  a   tungsten    substrate   was    sputtered   rather    than   the 
tungsten   itself.       The   results    of    this    experiment,    that   the 
sputtering   ratio,    atoms    removed   per    incident   ion,    increased 
with   increased   ion    mass    and   increased   ion    energy,    qualitatively 
suggested   an   ejection    mechanism.      A  few    years    later,    Von 
Hippel    and    Blechschmidtv     '    proposed    a    theory    which    described 
sputtering   as    an   evaporation   of    surface   atoms,    an    emission 
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mechanism.      Earlier,    Von   Hippel(3)    had  found  by   spectroscopic 
means   that   at  least   some   sputtered  atoms    were   in   an   excited 
state.      The   sputtering  theory   showed   that  atoms  in   the   region 
of   impact   could  rapidly   acquire   thermal   energy  if    the   kinetic 
energy  of   incident  ions   was    converted   to    thermal   energy  at   the 
target   surface.      If   an  atom   acquired  a   sufficient  amount  of 
thermal   energy,    it   would   then  evaporate  from   the   surface,    some 
atoms   evaporating  while  in  an   excited   state.      This   theory  was 
improved  by  Townes^'    in   194^*i    who   calculated  a   sputtering 
flux  and   the  number   of   atoms   evaporated  per   incident  ion. 
Keywell'5'    in   195kt    used  neutron  diffusion   theory   to   approxi- 
mate  atomic   interactions   within  the   crystal,    a  new   approach   to 
the   theoretical  investigation  of   sputtering.      Direct  application 
of   statistical   methods   to   sputtering  was   made  by  Harrison'"' 
who   envisioned   the   interaction   of   two   distribution  functions 
(the   crystal  lattice  and   the   ion  beam).      These   models,    based  on 
statistical   methods,    implicitly  accept   ejection   type   mechanisms. 

One  of   the   most   important   contributions   to   the   study   of 
sputtering   was   made  by  Wehner'''    in   1953  •      In   the  first 
sentence   of   his   paper,    Wehner   stated,    "The   most   widely 
accepted   sputtering  theory   is   the   evaporation  theory...",    but 
the   deposit  patterns   of   single   crystal  sputtering  which  he 
obtained   showed  pronounced,    high   density   areas   or    "spots". 
There   was   now   strong  evidence  for  a   momentum   transfer 
process,    and  further   development   of   the   evaporation   theory 
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ceased.       Shortly   after  Wehner's   findings    were   reported, 
Henschke^"'    proposed   a    theory   of    sputtering   based    solely   on 
classical    collision    theory,    treating   normal    and    oblique    incidence 
sputtering    separately.       Oblique    incidence    sputtering    could   be 
explained   by    two    body    collision    processes.       The    incident    ion 
penetrated    the    surface    layer    with    little    or    no    interaction,    was 
relected   outward   by    second   layer   atoms,    and    ejected   a    surface 
atom   by   striking   it    on  its    inside   hemisphere.       His    concept    of 
normal   incidence    sputtering   required    many-body    collisions    in 
which    the    ion    was    eventually    reflected    outward    to    sputter 
surface   atoms.       This    theory    was    plausible   for    oblique   incidence 
sputtering,    however,    the    case   for    normal    incidence    sputtering 
required    that    the   ion  be    reflected   inside    the    crystal.       (For    an 
ion   more    massive    than   the    target    this    requirement    cannot   be 
met.  ) 

Silsbee    *^'    noted    that    in   any   discussion   of    momentum    trans- 
fer   effects    in    sputtering,     the    geometry    of    the    crystal    should 
be    considered;    successive    collisions    in   a    crystalline    structure 
might   be    influenced  by    the    structure   itself.       His    calculations 
showed:    in   a   two— body    collision,    an   atom   has    a   departure   angle 
n.  with    respect    to    an    axis    specified   by    a    close    packed    direction. 
In    subsequent    collisions,    the    departure    angle  p^ ,    will   decrease 
if    the    energy   is    low    enough,    and    the    collision    sequence   occurs 
along   a    close    packed    row    of    atoms.       This    concept    of    "momentum 
focusing",     "Silsbee    chains",    or    "focusons"    was    considered   by 
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many   investigators,    excepting  Wehner^^),    to   be   a   fully    satis- 
factory   explanation    of    the    spots    in   deposit    patterns.       But    the 
deposit   pattern   is    only   one    of    many   observable   features    of 
single    crystal    sputtering. 

The    sputtering   ratio    is    another    important    characteristic 
of    the    sputtering   process.       If    a    theory    or    model   is    to    explain 
sputtering,    both    patterns    and    sputtering   ratios    must   be 

explained    qualitatively   and    quantitatively.      Almen   and   Bruce 

(  11 ) 
(AB)  measured    sputtering   ratios    of    a   variety   of    metals 

using   N ,    Ne,    Ar,    Kr,    and   Xe   ions    over   a    5-65   keV    range    of 
bombardment    energies.       (Although    polycrystalline    specimens 
were    sputtered,    the    qualitative    results    are    the    same    as    those 
obtained  for    single    crystal    sputtering.  )       They   noted    that    the 
sputtering   ratio    was    an   increasing   function   of  -  bombardment 
energy   provided    the   ion    mass    was   greater    than   the    target    mass. 
Sputtering   ratios    showing   a    maximum,    or    of    a    slowly   varying 
nature,    were    observed   in   all    cases    in   which    the    ion   was    lighter. 
Evidence    such   as    this    coincides    with   any   gross    conception   of    an 
ejection    mechanism. 

Single    crystal    copper    was    sputtered   by   argon   at   inter- 
mediate   energies   by   Magnuson   and    Carlson    ( MC )  *       '    and    South- 

(  13) 

ern,    Willis,    and   Robinson    (SWR)  .       MC    measured    sputtering 

ratios    (using    1-10    keV    argon)    for    the    (111),     (100)    and    (110) 
surfaces    and   found    that    sputtering   ratios    decreased   in    this 
order    of    surfaces.       SWR    sputtered    single    crystal    copper    with 
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1-5   keV   argon   and    measured   sputtering   ratios,    but   of    greater 

significance    was    the    quality   of    their    deposit   patterns    which 

clearly   showed   the   presence   of    only    certain   spots: 

(111)    surface  3-(110),       l-(lll)    but   not    seen   due    to 

beam    aperture. 
(110)    surface  1-(110),       2-(100) 

(  100)    surface  i+-(110),       1-(100) 

The   work   of   AB,    MC,    and    SWR    indicated   that    the    sputter- 
ing ratio    is   at   least   a   function   of    bombardment    energy    (or 
momentum)    and   of    mass    ratios.       However,    consideration   of 
spot   patterns   indicates    that   the    sputtering   process    itself    is 
also   highly   dependent    on   the    gross    and/or    surface    geometry   of 
the    crystal   and   on   ion   penetration   depth. 

Ion   penetration   of    solids   has   been   investigated    experiment- 
ally  and  by   computer    simulation.       Piercy,    McCargo,    Brown, 

(  lix  ) 

and    Davies      ^      investigating   channeling   of    various   heavy   ions 
in    monocrystalline   aluminum,    found   that    the   ion   penetration 
distance   increased   in   the    order    (111),     (100),    and    (110)    beam 
orientations.       This   is    the    same   order   as    the    sputtering   ratio 
decrease   found   by  MC.      A   number   of    studies    of    ion   penetra- 
tion have   also   been    made   using   computer    simulations    (Oen, 
Holmes,    and   Robinson'     5' ,    Robinson   and    Oen  ,    Harrison, 

Leeds,    and   Gay*  I '' )  • 

The    most   apparent    correlation   between   the   results    of    ion 
penetration   and    sputtering   studies    is    the    concept   of    trans- 
parency,   defined   by    Fluit,    Rok   and    Kistemaker'        '  .       Both 
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experimental   and    simulation   investigations    of    ion   penetration 
indicate    that    channels    exist   in    certain    crystallographic   directions 
of    a    crystal.       Results    of    normal   incidence    sputtering   experi- 
ments   indicates    that    sputtering   ratios    are   lower   for    surfaces 
in   which    these    channels    are   parallel    to    the   beam.       One    might 
initially    suspect    that    sputtering   ratios    will   be   low    if    ions    and 
target   atoms    are    confined    to    these    channels;    momentum    is 
directed   into    the   crystal   rather    than   laterally,    reversal   of 
momentum    to   produce    sputtering   not    occuring.      Accordingly, 
some   incident   ions    should   see   a   transparent   surface   rather 
than   a   wall   of    atoms. 

The   obvious    complexities    of    the    sputtering   process    may  be 
investigated   separately  by   digital   computer    simulation.       This 
of    course   involves    the    selection   of   a   gross    mechanism    (ejection 
or    emission)    and   a    model    (hardsphere   or    otherwise).      If 
factors    such   as    the   interatomic   potential   function   and   its 
parameters   are   known,    then   the   use   of    the   right    model   in   the 
computer    may  be    expected   to    simulate   the   actual    sputtering 
process.       Explanation   of    deposit   pattern   characteristics   and 
sputtering   ratios    will   necessarily  follow   from    an   exact 
simulation. 
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2.       PURPOSE    OF    THE    INVESTIGATION. 

The   belief    that    the   formation   of    spots    in    sputtering   deposit 
patterns   is    solely   a    consequence   of    momentum   focusing   was 
widely    accepted   for    a    number    of    years.       The    importance    of 
momentum    focusing,    not   only   in   spot   formation,    but   in   the 
entire    sputtering   process      is    now    believed    to    be    much    less    than 
originally   thought    liy_^->J  .       This   investigation   was    undertaken   to 
explore    the    single    crystal    sputtering   process   using   a   digital 
computer    to    simulate   a    copper    crystal    which   is   being  bombarded 
with    a   beam    of    argon   ions.       Primary    effort   was    directed 
towards: 

a.  Finding   the    mechanisms    which   result   in   atoms   being 
deposited    in    the    spots    or    high    density   regions    of    the 
deposit   pattern. 

b.  Determining   the    correlation   between   crystal    structure 
and    sputtering   ratios    and   patterns. 

c.  Investigating    the    surface   binding   energy   and   its    effect 
on   deposit   patterns. 
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3.       MODEL  c 

A.       A    space    lattice    is    established    whose    sites    represent 
the    equilibrium    positions    of    copper   atoms    in   a   face    centered 

cubic    crystal.      Atomic    spacing   is    that   determined   by   x-ray 

o 
crystallographic    studies    (for    copper,    a   =   3*  615   A).       The 

potential   function,     with    parameters   for    copper-copper    inter- 
action,   is    the    Born-Mayer    type    Gibson   Number    Two  .       The 
interaction   between   argon   and    copper    is    similarly    described 
using   a   Born-Mayer    potential,    V(r)    -    e           r,    but    the    numerical 
values    of    the    parameters   .    .    and   B  ,    for    this    function    are    those 
determined   by    Harrison,    Carlston,    and   Magnuson    (HCM)'2y) 
from    a    study   of    secondary    electron    emission.       The   range    of   both 
potential   functions    is    eroded   at   r0 ,    one   half    the   nearest   neigh- 
bor   distance    such    that    the   potential   and   force   go    to    zero    for 
atoms    whose    centers    are   farther   apart    than   2r    .       Use   of    this 
eroded   form    of    the   potential   functions    permits    establishment 
of    the    crystal   in   a    stable,    static    state    without   use    of    an 
attractive    potential.       This    equilibrium    state    is    destroyed    only 
when   an   atom    moves   from    its    lattice    site. 

The    lattice    as    a    whole    is    restricted   by    the    current    computer 
program    to    perfect    form.       Defects    such   as    interstitial   atoms 
cannot   be    used,    however,    irregularities    in    the    surface    layer    are 
introduced   by   using    two    types    of    variation  from    a    perfect 
surface;    vacancy   and    stub    conditions.       The   perfect    or    regular 
surface    is    identical    to   a   full   lattice   plane.       In   the    /acancy 
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configuration   the   atom    which    would   normally   be   hit   first  by 
the   ion   is    removed  from    the   perfect    surface.       The    stub 
surface   has   an   atom    added   on    the   perfect   surface   in   a    stable 
position   and   adjacent   to    the    target   atom.      Additional   atoms 
can   also   be   added   on   the    surface   or    removed   from    the   perfect 
surface    to    provide   a   variety   of    random    irregularities. 

B.       Each   ion   of    the   beam    is   approximated   by  a    single, 
neutral   argon   atom    whose   velocity   vector   intersects    the   plane 
of    surface   atoms    at   an  impact   point.       (The    term    ion  is   used 
throughout    this    thesis   for    the   incident   particle    to    avoid    con- 
fusion  with    target   atoms  .  ) 

Successive   runs   are    made   using   different   impact 
points    to    simulate   use    of    the    entire   beam    area.       The   use    of 
this    simple    model   for    the   ion  beam    is    subject   to   four 
conditions: 

1.  The   results    obtained   using   one   random    set   of    impact 
points    are    essentially    the    same    as    those   obtained 
using   any   other   random    set    of    impact   points. 

2.  The    majority   of   beam    ions    are   neutralized  before 
reaching    some   arbitrary  boundary   which    defines    the 
surface. 

3.  Prior    to    neutralization,    the    path   of    a   beam    ion   is 
not    appreciably    affected   by   the    surface   potential. 

4.  A    crystal   region   is   in   an    equilibrium    state    each   time 
a   neutralized   ion   impacts   in    that   region. 
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Ea  ch    of    these    conditions   has   been    satisfied    either 
in   the    model   itself    or   by    calculation.      A   full   discussion   is 
contained   in  Appendix  A. 
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4.       DYNAMICS. 

The   force   functions   for    copper-copper   and   argon-copper 
interactions    are   derived   from    the   respective   potential   functions. 
Atom    position   and   velocity   are    then   determined   using  Newton's 
Second  Law    of   Motion.       Normally,    one   would   determine    these 
values   by   numerical   integration    methods    over    small   time   inter- 
vals,   however,    this    process   uses    too    much    computer    time. 
Special    methods    approximating   integration  have   been   developed 
which   are   used   to    minimize    computer   run   time   yet    maintain   a 
good   approximation   to    the    integration. 

A.      Gay' 28)  developed   an   iterative    method   of    solution   which 
is    similar    to    that   used  by   Gibson,    Goland,    Milgram,    and 
Vineyard**"").       Both    Gay's    and    Gibson's    methods    use    the 
technique   of    replacing   differential    quantities    with   finite 
differences    to   approximate    the    equation    of    motion,    F/m    =    dv/dt 
\\\/  It   by   F/m    =  Av/^t.       Using   this    equation   and    the   relation- 
ships Ax   =   v/lt   and   v   =    (v   +   vQ)/2,    equations.  (3-1)    and    (3-2) 
are   obtained. 

Av   =   v   -   vQ        -         [f    (xQ)/m]    At  (3-D 

Ax   =   x    -   xQ         -       HF    (xQ)/m]    &  t/2   +   vi  At  (3-2) 

The   known   values    F(x    ),    v    ,    x    ,    and  At    in    these    equations 


o»       o 


are   usually   associated   with    a    time    t    =    t    ;     x   and    v,    the 

unknown  values,    are    to   be    evaluated   using    these    equations   at 

time    t    =   T(t    ,  £t):. 

Gibson's    technique   associates    the   value    of    x     with    a    time 

o 
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t        ™   t    .    but   v      is    associat.ed   with   a    time    t        -    t      -    A"t/2. 
ox  0  o  ov  o 

Current   values    of   x  and  v  are   alternately   computed   at.   intervals 
At   using   equations    (3-1)    and    (3-2). 

v(tQ    +  &t/2)    «   vQ  (tQ    -  &t/2)    +    ^F(xo)/ml(At 
x(t0    +tt)    =   x(tQ)    +J  [F(x0)/mW/2   +   v(tQ    +6t/2)Ut 
The   advantage   of    this    technique   over    that   of    evaluating  both 
variables   at   time    t   =■    t      +  &t   is    that   a    continuous    smoothing 
of    the    values    of   x   and   v   occurs    with    each   successive    computa- 
tion. 

Gay   recognized    the   need  for    some    sort   of    smoothing 
procedure    but    considered    that   an  averaging  process   for    equation 
parameters   rather     than   computed  results   would  result  in  a 
better   approximation.       His    method   replaces    the    evaluated 
force  function   F(x0)    in   these   two   equations   with   an  arbitrary 


force  function  f  ( x)  which   is   linear   in   some   interval   to  be 
determined.      Consider   the  following  two   step   cycle: 

Move   an  atom   from   x      to   x*    using   computed  values   of 
F(x0)    and   v0    to    solve   equations    (3-1)    and    (3-2).      Compute 
F(x^)    and   average    this    value   with    F(xQ).       Now,    move    the 
atom   from   Xj    to   xg   using   the   averaged  force,    F(xQ,  x.,  ) 
and  v0    to    solve    the    equations    using   the    same    constant   value 
of  At.      Let  f  now   be   defined  in   the   interval    (xQ,    Xo),    such 
that  f(x0)    "   F(x0'    xi  ^  *      Equations    (3-1)    and    ( 3-2 )    may 
now    be    written: 

v(t   +  At)    ■■■    v(t)    4     [i'(x0)/ml   A1 


x(t   +  Lt)    =    x(t)    +J[f (x0)/raj^t/2    +    v(t)(     &   t 
The   values    of   v   and   x   are   determinable   in   any   interval   in 
which   f    approximates   _F    to    some    arbitrary   degree,    therefore, 
in   every   complete    cycle    or    timestep   a   new    interval   and   new 
function   f    are    defined. 

A   determination   of    which    of    these    two    methods    is   better 
must  be   based   on   some    standard   or    condition.       In    many   cases 
the   accuracy   of    numerical    solutions    such   as    these    may   be 
compared    to    the    exact   value    one   obtains    by   integration.      But, 
there   are   no    truly    exact   results   in   a    simulation    since   a   physical 
model   is   in   the    computer   rather    than    the    computer   being  used 
as   a   means    of    solving   equations.       However,    if    the    model   is 
physical,    physical   laws    must   be    satisfied,    and   the    energy   gain 

or   loss    due    to    the    mathematical   approximations    must  be    small. 

(  31) 
This    condition   of    energy  balance   has   been   used   by   Harrisonv^    ' 

who   found    that,    in   simulations    involving   atomic   interactions    such 

as    channeling   or    sputtering,    use    of    Gay's    method   resulted   in   a 

better    energy   baiance. 

B.     The   value    of    these   approximation    techniques    may  be 

substantially   reduced   by   injudicious    selection   of    a   numerical 

value   for  &t.       Too    large   a   value   invalidates   the    approximation, 

but   too    small   a   value   increases    computer   running   time.       The 

program    used   in   this    simulation   of    sputtering   incorporates    not 

only   Gay's    approximation    method   but   also    a    procedure   which   he 

developed   for    automatic    adjustment    of  At.       Its   value   is 


calculated   using   the   relationship  6  t   -      or/v   "where  &t   is   replaced 
here  by   dt   to   indicate    that   its    value   varies   with   each   timestep. 
The   velocity,    v   is    the    magnitude   of    the   velocity   of    the   most, 
energetic   atom   in   that   timestep,    and    0  r    is   used   here   as   an 
input   parameter,    the    timestep    multiplier.       This    parameter   is 
assigned   the    computer   variable   name    DTI    and   is   defined:  the 
maximum    displacement   of    any   atom   in   any    timestep    shall  be    the 
numerical   value   of    DTI    in   appropriate   units   of    length.       Gay 
considered   the   value   of  &t,    consequently   that   of    DTI,    to   be   a 
function   of    ion   energy,    impact   point,    and    start  point,    but 
he   was   unable    to   find    the   relationship.       Trial   and   error    methods 
were   used   to  find   optimum   values,    and  reasonably  good  results 
were   obtained   using   these   values.      Johnson'  ^9 )    reported    that 
the   energy    balance   was    maintained   to    within  J>%   for   all  his   runs. 

The   difficulties    in    choosing   a   value  for    DTI   have   been 
largely   eliminated   in   this    study.       Some   problem    areas    such   as 
start  point  and  impact  point  dependence  were  avoided   by  always 
positioning   the   ion   tangent    to    the   first   atom    it   would  hit. 
This,    at   worst,    will   result   in  a    constant   error   for    which 
adjustments   may    be   made    ( this   has   not  been  necessary)  .      The 
positioning  procedure   is   described  in  Appendix  B  .      The   depend- 
ence on  ion  energy  is   inherent    to   the   DTI    computation  process 
and  Is    easily   explained   once    the    clutter    of    impact   point;:   and 
start  point,  dependence  has   been   cleared   away.       Use   of    the 
DTI    process    excludes   all   variable   dependence   except    the   Inverse 
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proportionality   to    velocity,    but   the   interval   in   which   the   force 
function   is    considered   linear    is    exactly   the    same    width   as    the 
numerical   value   of    DTI.      Accordingly,    a    straightforward 
analysis   is    made    to    show    the   interdependence. 

The   force   function,    F(r),    is    expanded   about   an   arbitrary 
point   r*   using   a   Taylor's    series   for   a   function   of    one   variable: 


F(r)=F(r*)    +  ~^F 

Tr" 


2 


(r-r*)    *    1  "^    F 


r# 


2!3r2 


r* 


(r-r*)     +...1     j  F 

k!lr^ 


(r-r*) 


Since    F(r)    is    derived  from    the   potential   function  V(r)    =    eA+Br, 

e  ■  -■-  r  ""S^  k 

3  ,     F(r)    is    necessarily   of    exponential   form    and  O    F     =   B    F 

B_  F-i    and    the    series    is    rewritten  in    the   following  form: 

F(r)-|F(r*)[l  +  B(r-r*i] UUr*)     1  2_    Bn   (r-r*)n    I   (3-3) 

The    quantity    (r-r*)    is   now    defined   as    a.  variable   displacement 
Or,    and    equation    (3-3)    may  be    written:    F(r)    =   f  (  C  r)    +   g(  o  r) 
where   f  (  0  r)    and   g(  o  r)    are    the    terms    in  braces    in   the   force 
function   expansion.       The   function   f (o  r)    is   assumed    to   be    the 
linear   function   used   in   Gay's    approximation    method;    the   function 
f(  o  r)    is    to    approximate   F(r).       Thus,    it   is   required   that    the 
ratio-\g(  a  r)/f  (  0  r)    be    some   fraction   less    than   one. 

Let   an   assumption   be    made   that   the    ratio    is    much   less 
than   one,    but    that    the   actual   value    is    the   fractional   deviation' 
from    linearity,    a   variable    whose   value   is    to   be    specified.       Then, 
since   B    is   a   negative    constant    (from    the   potential   function), 
the    truncation   error   for    the    series    is    easily  found,    and    the 


11 
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quantity   B  O  r     is    known   for    every    specified   value    of    fractional 
deviation.       But   B    is    constant,    and   once    the   fractional   devia- 
tion  is    specified,   Qjc    is   fixed.       Consequently,    specified    to    the 
fractional   deviation   from    linearity,    there   is   a   function  f(r) 
which   is   a   valid   approximation   to    the   function   F(r)    in   every 
interval    (r,    r+Or)    where   or    is   fixed.       The    quantity   6  r    is 
thus    the   value    to   be   used   for    DTI,    and  dt   =   dr/v   =    DTl/v. 

It   is    now    assurred   that  &t   will    always   be    the    maximum 
possible   value   which    minimizes    computer   running   time   without 
introducing   errors    due    to    the   non-linearity   of    the    true   force 
function.       One   would    expect    therefore,    that    the    energy 
balance    would   be    maintained    to    the    same    order    as    the   linearity 
of    the   force   function.       Further    consideration   will    show    that 
the    energy  balance    is    maintained    to    a    considerably   better 
degree    than    that   anticipated.       First,    since    DTI    is    the    maximum 
displacement    of    the    most    energetic   atom,    essentially   all   other 
atoms    will   be    displaced   a   distance   less    than   DTI.       Thus,    the 
interval   in   which    the   force   function   f(r)    is    used   is    smaller, 
and   f    is    an   even   better    approximation   to    F_.       Second,    for    a 
given   fractional   deviation   from    linearity,    the   value    of  or    is 
total   displacement,    not   x,    y_,    or    z    component   displacement j 
and   the   velocity   used    in    the    determination  fat.   -    DTl/'v   for 
each    timestep   is    the    total   velocity   v,    not    the    component 
velocity    vj  .       The   result   of    these    considerations    may    be 
summarized   by   a    calculation   for    the    maximum    displacement    oi 
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an   atom    along   the    i     component   of    its    displacement   vector: 
&Xj    =   v.&t,    but   v.  4    v  C   vmaxand  &t   -    DTl/vmax<       It   follows 
that  Axj    =    v;       (DTI,/vmax)   £    DTI. 

Tnese   analyses    appear    to    provide    a    tidy    solution   to    the 
time    optimization   and    energy   balance    problem.       This    would   be 
true    except   for    the    complications    introduced  by   the    eroded 
form    of    the   potential   function.       Consider    a    situation   in   which 
one   atom    is    separated   from    another   by   a   distance    ( 2r      +    d) 
where   d   is    some   distance   less    than   the    total   distance    D    that 
the   atom    will    move   in   the   next    timestep.      When   the    atom    is 
moved    the    distance    Dj  as    shown   in   figure    1,    the   force   on   the 
atom    at   its    new    position   is    the    same   force   it   would  have   if 
it  had    moved    only    (D-d)    units.       Since    the   force   is    zero   for 
r>£r0,    t^le    m°del   assumes    that   the   atom   had   a   velocity   v   = 
(D-d)/&t    (which   is    smaller    than   its    actual   velocity   v   =    D/&  t    , 
and   a    smaller   kinetic    energy   is    computed   using   this    smaller 
velocity.       This    difficulty    was    originally   foreseen   by    Gay   who 
included    corrections    to    the   force    calculations   for    such    situa- 
tions.      Both   Levy^30)    and   Johnson^"'    improved   the    methods 
of    corrections,    and   recent   improvements    in   the   force    calcula- 
tions  have   been    made   by   Harrison^    '  . 

The    continuous    improvement   of    the   dynamics    section   of 
the   program   has    resulted   in   a    model   in   which   not    only  has    the 
computer    running    time   been    minimized,    but   of    greatest   signifi- 
cance,   for    nearly   all    of    the    700-800    combinations    of    surfaces, 
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ion   energies,    and   impact   points    used   in   this    study,    the    energy- 
balance   was    maintained    to    within    1% . 
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5.       RESULTS. 

Reports    of    sputtering   studies    often   include   ad  hoc   formula- 
tions   of    mechanisms    which   are   used    to    explain   deposit   patterns, 
sputtering   ratios,    and    sputtered   atom    energy   distributions. 
The    sputtering   process    is    discussed    in    this    thesis    in    terms    of 
mechanisms    which  have   been   observed   in   the    simulation    to    sputter 
atoms .     The   observance    of    these    mechanisms    is    an   advantage 
peculiar    to   a    simulation.       Each    crystal   atom    (and   the   ion)    must 
be   identified  by   number,    at   least   for    purposes    of    computer 
calculations.       Since    every   atom    is   identified,    its    complete    track 
can  be   plotted   and   labeled.       The    tracks    of    selected   atoms    can 
then   be    superimposed    to    show    a    complete    set    of    interactions. 
The    mechanisms    which   were    observed    to    cause   an   atom    to 
sputter    are    considered   the   prime   observable    quantity   of    the 
simulation. 

The   results    of    the    simulation   are    discussed   for    each    crystal 
surface.       The    determination    that   an   atom   having   energy  _E   is 
sputtered,    was    made   using   a   probability-of -sputtering  function, 
P(E)    =    1    -    e"^/Eb    Where    E^  is    the    assumed  binding   energy. 
The   value   of    P  is    compared   with   a   random   number    R_  having   a 
value   between   0    and    1.       If    P  y  R ,    the   atom    is    considered 
sputtered.       (This    method   of    selection   has   been    compared   with 
one   using   a^step   function   in   which   all   atoms    with    E  >  E^    are 

assumed    sputtered.      A    selection    method   using   another   probabil- 

—  (  E  — E     ) 
ity-of-sputtering  function,    P(  E)    =    1    -    e  ~"    b      has    also   been 
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tested.      No  differences   were  observed  is  the   sputtering 
patterns  based  on  the  atom   selection   method.) 

Mechanisms   observed   to   cause  sputtering  are  discussed 
first  for   each  surface.      Deposit:  patterns,    sputtering  ratios, 
and  sputtered  atom   energy  distributions   are  then  discussed. 
The  deposit  patterns   shown  are   superpositions   of   those 
obtained  by  sputtering  the  crystal  separately  with  regular, 
vacancy,    and  stub   surface   conditions.      This  practice  follows 
that  used  by  ^Johnson'  ^'    who  found   that  pattern  features 
were  difficult  to  recognise  when  viewing  separate  patterns 
from  each  surface  condition.      The  energy  distributions  are 
also  superpositions  of  the   three  surfaee  conditions.      The  use 
of  superposition  is   considered  to   most  realistically  simulate 
the  condition  of  the  crystal  surfaee  at  various   times  during 
sputtering.      There  are  strong  arguments  however,    for  using 
only  the  regular  surface  condition.      Sputtering  ratios  are 
not  averaged  ( which  would  correspond  ft©  the  superposition 
©f  other  data)    since  a  probability  factor  for   the  existence 
©f  each  surfaee  condition  would  need  to  be  determined,      in- 
stead,  the  simulation  sputtering  ratios  for  vaeaney  and  stub 
surfaee  e©nditions  are  included  and  discussed  briefly  only  i 
completeness,      (The  (110)   surfaee  was  sputtered  only  at 
and  3  keV  using  all  three  surfaee  conditions.) 

C  Jlographie  nomenclature  is  used  gen<  rt   >|iy  in    t  n« 

diseussi©n  of  sp©es,     Reference  to  an  (hkl)   spot  does  ae1 


imply    (hkl/  ejection.       The^hkl^  directions    are   with   respect    to 
a    right  handed    coordinate    system    in   which    the    crystal   is 
described.       The    (hkl)    surface   is    contained   in    the    x-z    plane; 
the   y_   direction   is   into    the   drystal.       The    three    crystals    used 
in    the    simulation   are    shown    in   figures    2-i+;    all    contain    150 
atoms.       This    crystal    size   has   been  found   sufficient    to    contain 
nearly   all    energetic    collisions    for   ion   energies   up    to    7   keV. 

A.      The    (111)    surface. 

The    (111)    surface   was    sputtered   normally   with   argon   at 
1-5,    7,    10,    20,    and   1+0    keV.       (Although    the   potential  function 
is    considered   valid   only   for   ion   energies    less    than   about    7   keV, 
the   higher    energy   runs    were   used    to    search   for   additional 
mechanisms.)       Four    mechanisms    were   found.       These   were 
especially   evident    at    the   lower    (1-5    keV )    ion   energies.       Three 
of    these   are    classified    surface    mechanisms    since    only    surface 
atoms    were   found   to    participate   in   the   sputtering   event.       The 
fourth    mechanisms    is    a   deep    mechanism    in   which   an   inward 
moving   atom   is    reflected   and    then   initiates   a    sputtering    event. 
In   all    cases,    the    dominant    mechanisms    were   found    to   be 
surface    mechanisms. 

Except   for   head-on    collisions,    the   ion   is    scattered  by    the 
target    (atom    2)    with   a    component   of    momentum    parallel   to 
the    surface   as    well   as    a   much   larger   perpendicular    component 
into    the    crystal.       The    effect    of    this    combination   of    momenta 
is    to    drive   a   nearest   neighbor    (n.    n.  )  ,    for    example   atom    6    in 
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figure    5,    into    an   apparent  (112/    channel   formed   by   atoms    7, 
8,    and   9«       The    channel   is    apparent    since   it   terminates 
abruptly   thus    causing   atom    6   to    drive   atoms    7    and   8   into    the 
crystal   and   atom    9    outward.       Since   atom    6   passes    nearly 
directly  behind   atom    9>    the   impulse   is    more    normal    to    the 
surface    than   parallel    to    it.       This    is    termed   a    mole    mechanism 
since   one   atom   burrows   between    two    layers    parallel    to    the 
surface    to    sputter    an   atom    in    the    outermost   layer.      Atoms 
sputtered   by   this    mechanism    are   almost   always    a   next   nearest 
neighbor    (n.n.n.)    to    the    target   and   located   in   or    near    the 
sextant   defined   by   the   impact   area.       They   usually   sputter 
with    greater    than    10    ev.       The    target   always    receives    the 
majority   of    energy    transferred   by   the   ion,    but   it   is   normally 
driven   into    the    (111)    trigonal   array   of    atoms    directly  behind 
it.       These    three   atoms    act   as    buffers    and   dissipate    the 
target's    energy   and    momentum    into    the    crystal.       The    target 
is   reflected,    but   it   does    not   retain   sufficient    energy    to 
sputter   a    surface   atom    nor    sputter    itself. 

The    second   and    third    surface    mechanisms    occur    when    the 
impact   parameter    is    about    one    third    of    an   atom    radius    or 
less.       In  both    mechanisms    a   n.n.,    atom   6,    is    struck   by    the 
ion    such    that   the   parallel   and   perpendicular    momentum    compon- 
ent   magnitudes    which    it    acquires    are   of    the    same    order.       ii 
the   perpendicular    component.:   is    greater,    atom    8    scoops    a   n.n.n., 
atom    9,    relatively   high    up   in    the    surface    plane    of    atoms    and 
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into    their    outer   hemispheres.       These   atoms    are   densely   packed 
in    the    surface    and    act   as    a    rigid   reflector    for    atom    9.       It    is 
found    that    atoms    sputtered   by    the    scoop    mechanism    have 
ejection   angles    nearer    the   normal    than   would  be   found   in    the 
absence    of    reflection    in    the    surface.       If    the    perpendicular 
momentum    component    of    atom    6    is    nearly    the    same    or    less 
than    the    parallel    component,    it    squeezes    atom    9    against   its 
surface   neighbors.       The    squeezing    causes    the    surface   plane    to 
warp   and   atom    9,    10,    or    11    is    sputtered.      Atoms    sputtered 
by    either    the    scoop    or    squeeze    mechanism    usually    sputter    with 
less    than    10    ev.      Atoms    sputtered  by   combinations    of    the 
three    surface    mechanisms    are    found    to    sputter    with   higher 
energies    than   when    sputtered  by   a    single    mechanism. 

The   fourth    mechanism    requires    that    an   atom   be    reflected. 
The    atom    which    is    most    often    apt    to   be    reflected   is    the 
target.       If    the    impact    parameter    is    greater    than    one   half 
the    atom    radius,    the    target    is    driven    towards    an    edge    of    the 
trigonal    array    behind    it    rather    than    into    it.       It    will    penetrate 
this    edge    if    it   has    sufficient    energy,    but   it    will    lose    most    of    the 
energy    during    the    penetration.       The    reflection    occurs   from 
third    layer    atoms    and    the    target   atom's    reversed    momentum    is 
transferred    to    a    second    layer    atom.       The    second    layer    atom 
will    then    sputter    one    or  two    surface   atoms    at    energies    upwards 
of    5    ev.       At    lower    ion    energies    the    target    atom    cannot 
penetrate    the    array    edge    and    is    reflected   from    the    second 
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layer;    it    may   or    may    nor    sputter   a    surface   atom,    but   it 
rarely    sputters    itself. 

When   the    crystal   -was    sputtered   at    7   keV ,    the    ion   was 
found   to   penetrate    well   into    the    second   layer.       This    is   a 
significant   increase   in   ion   penetration   depth    since    5   keV    ions 
are   found    to   penetrate    only    to    the    order    of   half   an   atom    radius. 
The    effect   of    deeper   ion   penetration   is    to    cause   a   lateral 
compression   of    the    second   layer    with    accompanying   warping. 
Second   layer    atoms    which   are    squeezed   outward   will    sputter    a 
surface   atom   but   are   rarely   found    to    sputter    themselves.      At 
the   higher    bombardment    energies    tested    (10,    20,    l+®    keV),    the 
number   of    second   layer    atoms    which    sputtered   increased   nearly 
proportionally   to    the    increase   in   ion   energy.       The    sputtering 
process   at    these   higher    energies    appears    to   be   predominantly 
by   the   three    surface    mechanisms    described  for   lower    energy 
sputtering.      Mechanisms    peculiar    to   high    energy   sputtering  have 
not  been   observed   although    this    may   be   due   to    the   limited 
energy   range   for    model   validity. 

All    (111)    simulation   deposit   patters    showed    the    character- 
istic  features    of    three    (110)    spots    and   a    central    (ill)    spot. 
The   2   keV   pattern   definitely   showed    the   presence   of    three 
additional   spots    which    were   found   in    the   regions    in    which 
streaking   occurred   at   other    energies.       These   patterns    are 
shown   in  figures   6-21.       The   appearance   of   a  hexagonal   pattern 
at   2   keV    was   a    surprising  result,    but   it   is   not   a   unique 
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occurrence   in    sputtering    studies.      Anderson   and  Wehner         ' 
found   a  hexagonal   pattern   for    (111)    copper    sputtered  by 
mercury.       The    expected    trigonal    pattern    was    found   for 
energies    up    to   2+00    ev,    a   hexagonal   pattern   appeared   at   2+00 
ev,    and   it   disappeared   as    the   ion   energy   increased.      Very 
recently,    Robinson   and    Southern^    '    have   found   additional 
spots    near  ^  114^  positions   for    (111)    gold    sputtered   "with 
2+   keV    argon. 

The    pattern   features    vary    slowly   with   binding   energy 
provided   its   value   is    restricted    to    the   range   2. 50-3*50    ev. 
The    upper    limit    is    the    sublimation    energy   for    copper    (33)> 
above   it    the   patterns    appear    to    deteriorate.       This    was    most 
evident   at   a   3    keV   bombardment    energy.       A   set   of    patterns 
at   this   bombardment    energy   is    shown    in   figures    8-12+   for    a 
binding   energy   range    of    1. 50— 5*50    ev.       The   pattern   deteriora- 
tion  is    most    evident    in    the    region    of    the    (111)    central    spot. 
Figures    11    and    12    show    the    transition   from    a    well    defined 
central    spot   at   3*00    ev   to    complete    deterioration   of    the 
central    spot   at   2+.00    ev. 

The   half -intensity    width    of    the    (110)    spots    is    estimated 
from    pattern   data   and   from   numerical   data    to   be    about    11 
degrees.       The    width    of    the    central    spot    is    an   unreliable 
datum    and    was    not    determined. 

Sputtered   atoms    were   always   found    to   be    surface   layer 
atoms    for    sputtering   at    5    keV    or    less.       The    atoms    which 
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were    sputtered    most   frequently    appeared    to    naturally   group    by 
energies    into    two    categories;    atoms    with    perpendicular    energy 
greater    than   3    but   less    than    10    ev   and   atoms    with    perpendicular 
energy  greater    than    10    ev.       The   higher    energy   atoms    were 
found    to   be    equally   distributed   between   the    three    (110)    spots 
and   the    central    spot.       They   were   also   found   in    the   narrow 
sectors    defining    the    streaks    and   at   distances    from    pattern 
center    corresponding   to    the   distances    of    the    (110)    spots. 
This   distance   was    found    to   be    0.7    units,    the   distance    measured 
for    the    spots    in    the   patterns    of    SWR.       The   lower    energy 
atoms    were   usually  found   in    the    central    spot    which    explained 
its    sensitivity   to   binding   energy. 

The    correlation   of    sputtering  frequency   and    crystal 
location    of    sputtered   atoms    was   found   using  frequency-location 
diagrams.       These    are    shown   in   figures    22-29   for*    regular 
surface    sputtering.       Similar    diagrams   have   been    made   for 
sputtering   the   vacancy   and    stub    surfaces.       It    was   found    that 
a    1    keV    ion   does    not    cause   frequent    sputtering   of    n.n.n:. 's    or 
n.n.'s    in    the    sextant    defined  by   the    impact   area.       However, 
at   2    and   3    keV,    n. n.n.'s   are   found    to   be    sputtered    most 
frequently,    and   at   higher    energies,    both    n.  n.n.'s    and   n.n.'s 
are   frequently    sputtered.       The    three-fold   relationship   between 
the   frequency   of    sputtering,    location    of    the   atom    with    respect 
to    the    target,    and    the   ion    energy    may   be    clarified   by    consider- 
ing  the    sputtering    mechanisms. 
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At    1    keV    and    an    impact    parameter    of    about    three   fourths 
an   atom    radius,     the    ion    is    found    to    penetrate    less    than   half 
an  atom    radius.      At    energies    up    to    5    keV    the   penetration   is 
only    slightly   greater.       It    was    found    that,    although    the    target 
is    always    driven   into    the    crystal,    the    n.n.    which    initiates    the 
majority    of    sputtering    may    not    be    driven   far    enough    into    the 
crystal    to    provide    the    scoop    mechanism.       If    it    is    not    driven 
inward   a    sufficient    distance    to    scoop    the   n.n.n.,    its    energy 
will   be    propagated    in    the    surface    along    close    packed    rows 
which    originate    at    the    n.n.n.       This    is    not    a    "focuson."      At    an 
early   point    in    the    propagation    the    surface    will   have    warped 
sufficiently    to    cause    an   atom    to    sputter.       The    sputtering   is 

a    result    of    a    squeeze    mechanism    but   not    one    directly    involving 

I 

the    n.n.       As    the    ion    energy   is    increased    the    scoop    and    rnole 

mechanisms    are    more    apt    to    occur    although    the    squeeze 


mechanism     is    found    to    occur    at   all    ion    energies.       The    transi- 
tions   from    dominance    of    one    mechanism    to    another    with 
increasing   ion    energies    may   be    inferred    from    the    profiles    shown 
in   figure    31    for    atoms    which    sputter    most   frequently. 

A    relatively    constant    sputtering   frequency    is    observed   for 
atoms    102    and    132,    but    the   percentage    of   high   and   low    energy 
sputs    varies    for    atom    102    which    is    a    n.n.n.       The    percentage 
does    not    vary    for    atom    132    which    is    neither    a    n.n.n.    nor    a 
n.n.       Atom    102    sputters    most    frequently    at    low    ion    energies 
by    the    squeeze    mechanisms,    at   higher'    ion    energies    by    the    scoop 
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and    mole    mechanisms.       These    relationships   have   been    confirmed 
by   the   atom    track   displays.       They   also   have    shown   that   atom 
132    is    sputtered   by   atom    102    through   a    squeeze    mechanism 
alone   at   low    ion   energies   but   in    conjunction   with    the    mole 
mechanism    at  higher   ion   energies.       The   profile    shown   for    atom 
86   is    the   result   of    simple   reflection   from    second   layer   atoms. 
Atoms   87    and    101,    which   are   not   in    the   impact    sextant   but 
are    sputtered   frequently,    are    sputtered   by   the    squeeze 
mechanism . 

Figure   32    shows    profiles    of    atoms    sputtered  by   the   deep 
mechanism.      At   low    ion   energies    the    target   reflects   at   an 
oblique    angle   from    second   layer    atoms    rather    than   penetrating 
the   layer,    and   it    enters    one    of    the   apparent  ^112^   channels. 
Atoms    25    and/or    55    are    then    sputtered   directly   by    the    target. 
(Although    this   is    similar    to   a    mole  mechanism    it    is    a   distinct 
mechanism    since    the    target    must  be   reflected   in    this    case   but 
not   in    the   former    case.)      At   higher    ion   energies    the    target 
will   penetrate    the    second   layer    and    is    reflected   from    third 
layer    atoms.       It    does    not    channel,    but   it    causes    a    second 
layer    atom    to    sputter    surface   atoms.      Accordingly,    atom    55    is 
sputtered   by    the   deep    mechanism    more    often   than   atom    25    is 
sputtered   by   a    channeled    target.      Atom   40      or    70      is 
sputtered   by   a    squeeze    mechanism,    usually   in   conjunction   with 
the    sputtering   of    atoms    25    and    55* 

These   results    for    the    sputtering   of    a   regular    surface 
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have   been   found    to    be    generally   applicable    to    sputtering   vacancy 
and    stub    surfaces.       Equal   numbers    of    atoms    from    the    three 
surfaces    are   found    in    the    spots   and   in    the    streak   regions. 
It   is    also   found    that    these   atoms    are   usually    the    same   from 
all   three    surface    conditions.       They   differ    only   in   energy   since 
it   is    the   ion   penetration   distance   which   indirectly   determines 
the    magnitude   and    direction   of    momentum    transferred    to 
surface    atoms.       The    vacancy    surface    sees    an    ion   -which 
penetrates    deeper;    this    is   analogous    to    an   ion   of    greater    energy, 
In   the    stub    condition,    the   ion   does    not   penetrate   the    surface; 
it    transfers    energy   to    the    stub    atom.       The    stub    then   assumes 
the   role    of    the   ion   impinging   on    the    surface,    but    the    stub 
appears    as    an    ion    with    less    energy. 

Figure   33    shows    the    sputtering   ratios   for    each    surface 
condition    as    a   function    of    ion    energy;    assumed   binding    energy 
is    a   parameter.      A   statistical   variation    of    +  20$    of    the 
sputtering   ratio   has   been  assumed.       This    would   normally  be    an 
unreasonably   large    deviation  for   laboratory   results,    but   it   is 
considered    conservative   for    the    numerical   results   from   a 
simulation   in   which    many   parameters    are   unknown.       The   actual 
value    of    sputtering   ratios    obtained   for    the   regular    surface    is 
generally   low    by   2-3    atoms/ion   for    a   binding    energy    of    3-5    ev. 
Further    investigation    showed    that,    in   a   few    cases,    surface 
atoms    3    planes   below    the    impact   point    would   be    sputtered   with 
about    4-5    ev.       The    average    increase    in    the    (ill)    sputtering 
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ratio    over    the    1-7    keV   range   is    estimated    to   be   no    greater 
than   2.0    atoms/ion  for    a  binding   energy   range    of    1.5    to   3«5 
ev.       This    range   is    shown  by    the    dotted   lines    in   figure   32a. 
Vacancy   and    stub    sputtering   ratios    show    erratic   behavior    over 
the   energy   range    studied. 

Energy  distributions    of    sputtered   atoms   are    shown   in 
figures   35-39.       The   preponderance   of    atoms   with   energies    less 
than   5    ev   is    caused  by   the   lack   of   an  intrinsic    surface  binding 
energy.      The   sputtering   selection  process    considers    all  atoms 
having  _E   greater    than   1    ev.       Thus    the   P(0.50)    energy,    using 
an   assumed  binding   energy   of   3«5    ev,    is   2.42    ev.       This   is 
especially   evident   for    the    1    and   5   keV   distributions.      If    this 
region   is   disregarded,    the   peaks    which   occur   at   5-8   ev  for 
all   ion   energies   are   assumed   to  be    the   maximum   for    each 
distribution.      Secondary    maxima   appear   at    16    ev  for    1   keV 
sputtering,    16    ev   for   2   keV ,    and   at   12+   ev  for   3   keV    sputter- 
ing.     At  higher    ion   energies,    a   secondary   maximum    may  be 
present   at   2+7    ev  for    5   keV    sputtering.       The   number   of    atoms 
sputtered   in   the    simulation   is    too    small    to    make   a   definite 
statement   regarding   maxima.       Differences   are   usually    measured 
by  one,    at   most,    two   atoms. 

B.     The    (100)    Surface. 

The    (100)    surface   was    sputtered   normally   with    1,    3»    5» 
and  7   keV   argon.      Mechanisms    observed   to   cause   sputtering 
were   found    to   be   nearly   identical   to    the    surface    mechanisms 
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discussed   for    (111)    sputtering.      A   mechanism    directly- 
comparable    to    the   deep  mechanism    was    not   observed.       Instead, 
it    was    found    that    the    scoop    and    mole    mechanisms    are    enhanced 
by    the    presence    of    the    (110)    channels    parallel    to    the    surface. 
Surface    atoms    which    are    driven   into    these    channels    do    not 
have    to   burrow   between    the   first   and    second   layer    in  order 
to    sputter   nearby    surface   atoms.      Additionally,    it   was   found 
that    these   atoms    will   frequently   sputter    themselves   by 
reflection   from    second    layer    atoms.       The    squeeze    mechanism 
was    observed    to   be    the    most   effective    sputtering    mechanism 
at   3    and    5   keV.       At    1    and    7    keV    ion   energies,    variations    of    the 
scoop   and    mole    mechanisms    were   dominant.       The   low-high   and 
middle   ion    energy   dependence   observed   for    the   dominant 
mechanism    is    directly   related    to    the    shape    of    the   impulse   by 
which   the    tion   transfers    energy   to    a    n.n. 

At    low    ion    energies,    the    impulse    is    sufficiently   broad    that 
the    n.n.    is    directed    into    the    edge    of    the    square    of    atoms 
behind    it.       The    n.n.    will    enter    the    (110)    channel    without 
immediately    sputtering   a    surface    atom.       When    the    impulse    is 
narrow,    the    peak   force    is    not    necessarily   greater.       The    target 
has    received    most    of    the    energy   given    up   by    the    ion,    and    the 
ion-n.n.    impact   parameter    is    larger.       The    n.n.    will    then   scoop 
or    squeeze    a    surface    atom    causing    it    to    sputter.       At   high    ion 
energies    (7    keV)    the    ion   was   found    to    penetrate   at   least   into 
the    second   layer    where   it   was    deflected    towards   but   not   into 
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a    (110)    channel.      This   was    observed   to    cause    second   layer    atoms 
to   be   scooped  up  and   sputter   surface   atoms.      Second   layer   atoms 
are   also    squeezed   by   their   neighbors   and    may   sputter    through 
the   vacancy  left  by  a   sputtered   surface  atom. 

Sputtering  at  3   and   5   keV   was   found   to  be   caused  by   the 
scattering  of   surface  atoms   along   the   surface   as    the  ion  was 
being  reflected  from    second   layer   atoms.      The   squeeze 
mechanism    was   observed    to   sequentially  sputter   atoms   along 
a   close  packed   surface  row.      The   occurrence   of   a   mole 
mechanism    was   conspicuously  rare;    even  at   5   keV   the   energy 
propagation   was    clearly  restricted   to   propagation  parallel   to 
the   surface   layer   of   atoms   and   to   propagation  into    the   crystal. 
The   only  definite   momentum   reversals   observed   were  for    the 
target  atom   and   the  ion. 

The   characteristic  features   of    (100)    sputtering  deposit 
patterns   are  four    (110)    spots   and  a   central   (100)    spot.      The 
patterns   are  usually   outlined  by   a  hypocycloid-shaped  hagej 
the   (110)    spots   form   the   cusps   of    the  hypocycloid.      Simula- 
tion patterns   were  found    to   show    these   characteristic 
features   quite  well  for   3   and   5   keV   sputtering.      The  3   keV 
pattern    (figure   i+1)    is   a  very   good   likeness   to   one   at:  eJ .  b   keV 
reported  by  SWR,      The    (1.10)    spot  distances  from   patti  m 
center   were  found    to  be   1.0    TC   unit   which   corresponds    to 
the  distances   In   the   2.5   keV    experimental  pattern.       Ii     was 
further  determined   from  numerical  data   that   the   (110)    spots 
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were   actually   elliptical   rather    than    circular.       This   latter 
characteristic   is    undoubtedly    the    result    of    distortion   inherent 
in   flat   plate    collection   of    sputtered   atoms    v    J '  .      The   3    keV 
pattern   was    found    to    remain   essentially   unchanged   as    the 
binding   energy    was    varied   from    2.50    to   3  •  50    ev.       The    5   keV 
pattern    (figure   42)    was    almost   identical    to    the    one   at   3    keV 
only   when   a   binding    energy   near    2.00    ev   was    used.      As    the 
binding   energy   was    increased,    the    (110)    spots   became   less 
well   defined    (figure   2+3)    and    spots    corresponding   to  ^  210^ 
became    the    most   prominent   feature.       The   intensification   of 
^210^  regions    was    accompanied   by   a   loss    in   definition   of    the 
hypocycloid   outline.       Similar    results    were    obtained   when    sputter- 
ing  at    7    keV.       The   hypocycloid   outline    in   the   pattern   could 
only  be    observed   for    a  binding   energy   near    1 .  50    ev. 

The   deposit   pattern  for    1    keV    sputtering    (figure   2+0 ) 
showed   very   little    similarity    to    those    obtained   from   higher 
energy   sputtering.       Four    spots    corresponding   to  ^211^  were 
observed;    four    (110)    spots    may   also   be   defined  but   they   are 
extremely   diffuse.       (The    (110)    spots    are   best    seen   in   the 
point   plot   in  figure   2+0a.)       The   finding   of    (211)    spots   at 
1    keV,    and    the   appearance   of    (210)    spots   for    5   and   7    keV 
ion   energies    suggested    that    the   pattern    might  be   rotated   45 
degrees   for   low    ion   energy   sputtering.      The   possibility   that   a 
rotation   occurs   has   been   investigated  t   but   one  has   never 

been   experimentally   observed. 

45 


The    different    (100)    pattern   features    which    appeared    with 
varying   ion    energies    and    binding    energies    were    a    marked    contrast 
to    the    relative    constancy    observed   for    (ill)    patterns.       (The 
use    of    a   binding    energy    of    3«00    ev    resulted    in    comparable    (111) 
patterns   for    all    ion    energies.)       The    lack    of    constancy   for    the 
pattern    simulations    may   be    explained    as    a    result    of    the    sub- 
traction  of    a   binding    energy   from    the    perpendicular    component 
of    a    sputtered    atom's    energy.       (This    is    done    to    simulate    the 
energy   lost   in      overcoming   the    surface    potential.)      An   atom 
sputtered    with    small    perpendicular    energy,    less    than    15    ev    for 
example,     may    suffer    an    apparently    small    change    in    its    perpendic- 
ular   velocity,    but    the    direction    of    its    velocity    may    be    consider- 
ably   altered.       The    average    perpendicular    energy    of    atoms    found 
in    spot    regions    is    about    15-20    ev   from    sputtering    the    (100) 
face.       It    is    about    20-25    ev   for    (111)    sputtering   for    which 
variations    of    up    to    1    ev   binding    energy   have   had    little    effect    on 
pattern   features.       This    explanation    does    not,    however,    explain 
why   the    pattern    at    3    keV    is    valid   for    a    range    of    binding 
energies.      A    consideration   of    the   atoms    sputtered    may   clarify- 
but   not    completely    explain    this    situation.       The   dissimilarities 
seen    in    the    (100)    patterns    may    be    directly    inferred    from    the 
frequency-location    diagrams    for    (100)    sputtering    (figures    46-49)- 

When  the  energy  propagation  is  restricted  to  the  surface, 
as  it  is  for  1  keV  sputtering,  the  (ill)  surface  can  dissipate 
energy    through    four    atoms    in    the    momentum    forward    semi- 
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circle.       The    (100)    has    only   three   atoms    comparably   located 
to   dissipate    the    energy.       Track   displays   have    shown   that   for 
1    keV    sputtering,    atoms    26,    116,    and    132    are    sputtered   by   a 
combination   of    mole   and    squeeze    mechanisms.       The    mole   is 
predominant   in   sputtering   the   atoms   nearest   the    target,    the 
squeeze    mechanism    becoming   predominant   for    surface   atoms 
■which   are   located   further   from    the    target.      Atom    26    sputters 
less   frequently   than   atom    116    since    it   is    on   the    opposite    side   6'c 
of    the    crystal   with    respect    to    the   impact   area.       Both   of 
these   atoms    are   found   in    (211)    spots,    but   atom    116    is   also 
found   in    (110)    regions    of    the   pattern.      Atom    132   is   found   in 
(211)    spots    and    the    (100)    spot;    it   is    sputtered   only   by    the 
squeeze    mechanisms.      When    the   bombardment    energy   is   3    keV, 
the   deposit   pattern  has    its    expected    characteristics.       Correla- 
tion  of   high   atom    density   regions    with   frequency-location   data 
showed    that    the    (110)    spots    are    the   result   of    sputtering   of 
n.n.n.'s.       Concurrently,    the   hypocycloid    outline    was   found    to 
be   formed  by:     (l)    the    sputtering   of    the    same    atoms    which 
formed    (211)    spots    at    1    keV ,    and    (2)    the    sputtering  of    the 
n.n.n.'s   in    the    surface.       These   atoms    were   not    densely 
deposited   in    the    (211)    regions    but   formed   diffuse    ellipses    whose 
semi— minor    axes    were   along   quadrant   bisectors.       The   appearance 
of    well    defined    (110)    spots    at   3    keV   is    strongly    suggested  by 
the   3    keV   frequency-location   diagram,    figure   2+7.      Atoms    70 
and    72    were    observed    to   be    sputtered   with   about    the    same 
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frequency  as   atomy    55,    5b   and   8b,      Sine©   atoms    ?0    and 
are   along  a    (100)    axis,    one   might   expect    bhai     fehey    would   be 
found  in   (110)    spot!,    and   this  has  been  observed.      Atoms   55 
and   56  are   in  positions   relative    be    the   target   such    that   they 
would   sputter'    within  adjacent  kb   degree   sectors   in   quadrants 
II   and   Hi.      These   atoms   are  found   in   the  hypoeyeloid   outline* 

The  absence   of   the  hypoeyeloid  at   1    keV   and   its   presence 
at  J,    3   and   7  keV    is  readily   apparent  when   one   considers   the 
sputtering  profiles   in  figure   50.      At   1   keV   ion   energy,    n.n.'s 
are  never   sputtered,    however,    they   are   sputtered   with 
increasing  frequency  as   the  ion  energy  increases.      Atom   86, 
the  n.n.    In  the  quadrant   containing  the  impact    area,    exhibits 

this   behavior      if   increasing  sputtering  frequency,      Ai ?g 

shows  a   sputtering  profile  which   indicates   the    (110)    6pe1 
formation  at  ion  energies  of  3  keV   or  greater.      The  constant 
nature  of  the  profile  for   atom   IQg  at  3   keV    and  higher   ion 
energies   is   indicativi        i     bhi     bypooyelold   outline   rather   than 
(211)  spot    formation*      Atoms   Vhf    116  and   1JI   shew   a   generally 
decreasing   frequeney  Of   sputtering   with  high  ion  energy.,      This 
is  eonsistem    with   bhs  loss  of  (§11)   spots  at  higher  Ion  energl 

The   braek  patterns  observed  for  5  and  7  keV    (100) 
sputtering  showed    bhai    atoms   ejected   by  a    squeeze   m    i  h  mil  R) 
are  often     |eoi    d  in  direotions  opposite  bo  bhese  one  would 
■  ■/,'/.  m-  solely   ■//(   i.-jh    basis  of   bhe  Leeatien  <<\    bh<    atnm   with 
respect  bo  the  Impaei    pointi     These*  ■oeourrenees  were  sai 


by   two   factors:     (l)    the    atom    was    squeezed   against   its 
neighbor    and    reflected    outward    and    with    its    parallel    component 
of    momentum    reversed    rather    than   being    strictly    ejected    out- 
ward.      (2)    Sputtering  by   near    simultaneous    squeeze   and    mole 
mechanisms    results    in    ejection    of    an    atom    where    the    ejection 
direction   is    dependent    only    on    the    impulse    delivered   by    the 
channeling    atom    in    the    mole    mechanism. 

Sputtering   ratios   for    the  (100  )'■  surf  ace;are  shoWnun  figure  51  • 
The    correspondence   between   simulation   values    and    experimental 
values    is    quite    good   for    the    regular    surface   for    a   binding 
energy    of    3*50    ev.       The   ratios   for    vacancy   and   stub    surfaces 
show    a    closer    correlation    to    the    regular    surface    for    (100) 
sputtering    then   for    (111)    sputtering. 

The  energy  distributions  of  sputtered  atoms  are  shown 
in  figures  52-55.  No  specific  maxima  other  than  that  near 
5-7    ev    appeared    to    be    present. 

C.      The    (110)    Surface. 

This    was    the    last   of    the    three   face    centered    cubic    crystal 
surfaces    to    be    sputtered    in    this    simulation.       It    is    purposely 
the   last    to    be    discussed.       First,    experimental   deposit   patterns 
from    (110)    sputtering    show    only    a    large    central    oval   area    for 
bombardment    energies    greater    than    a    few    hundred    ev    (  10 )  . 
Second,    a    unique    feature    of    this    surface    is    the    (110)    channel; 
no    other    (hkl)    surface    in    a   face    centered    cubic    crystal    shows 
^hkl^    channels.       The    sputtering    mechanisms    occurring   for    1,    3» 
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5    and   7    keV    ion   energies    were    most   readily   observed   for 
sputtering   of    this    surface,    and    they   are    interrelated   with    the 
presence   of    (110)    channels.      The    mechanisms    are,    again, 
identical   in   concept    to    those   previously   discussed   for    (111) 
and    (100)    sputtering.      Before    discussing    these    mechanisms 
separately   in   terms    of   individual    crystal   atoms,    the   features 
of    the   frequency-location   diagrams    (figures    56-59)    are 
summarized. 

Sputtering   of    the    (110)    surface   differs    significantly   from 
the   sputtering   of    (100)    and    (111)    surfaces:    (1)    the    target 
atom   was    the   atom    sputtered   most  frequently  at  all   ion 
energies.       (This   is    in   agreement   with   results    obtained   by 
Levy  .)       ('d)    the   number   of    sputtered   atoms    with   high 

energies   were   generally   those   in    the    close   packed   row    contain- 
ing  the   target  atom,    but    with   the.  target  atom   located  at 
the    center    of    the   row    rather    than   at   the   origin   of    the   row . 
These   findings    are    wholly    consistent   with    expected   results 
when  one   considers    the   sputtering   mechanisms   with  respect-   to 
this   particular   surface. 

The  frequent   sputtering  of   the   target:   is    made   possible   by 
the   nature   of    the   squeeze  mechanism    but    the    mechanism    does 
not   itself    cause   the    sputtering.       When   the    target:   is   driven 
into    the    crystal   It    squeezes    the   n.n.    (in   its   row)    and   ejects 
it.      Once   the   n.n.    position  is   vacant,    the   target    will    travel 
through   a   large,    potential-free   area   before    striking  a   second 
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layer    atom.       The    target    strikes    this    second   layer    atom    and   is 
free    to    reflect   outward   without  finding   a    surface   atom    directly 
in   its   path.       The    target    then    transfers    most   of    its    parallel 
momentum    to    a    n.n.n.    in    the    row.       The    parallel    impulse    is 
then   propagated   down    the    surface   row.      Warping   of    this    one 
row    causes   additional   atoms    to    sputter.       This    sputtering 
sequence   is    also    observed   as    a   result   of    the   ion  striking   the 
n.n.    which   is   in   the    quadrant   defined   by   the   impact   area.      A 
rather    surprising   continuation   of    this    sputtering   mechanism 
was    observed    to    occur    in   atom    rows    located  both   above   and 
below    the   horizontal   row    containing   the    target.       Second   layer 
atoms    which   receive    energy  from    the    target   or    ion   are   driven 
into    (111)    apparent    channels.       They   pass   behind   another    second 
layer   atom    causing   it   to   be    ejected   outward   and   strike    two 
adjacent    surface    atoms.       Sputtering   is,    again,    initiated   in   a 
surface   row.       Sputtering   profiles    for    the    target    (atom    72) 
and   atoms   2+2 ,    71,    73    and    102    are    shown   in   figure   60 .      Atoms 
42    and    102    are    the   n.n.'s    to    the    target;    atoms    71    and   73    are 
those   which   initiate    sputtering   in   the   horizontal   rows    containing 
them.       Second    layer    atoms    are    sputtered    more   frequently    as 
the    ion   energy   is    increased. 

Simulation   deposit   patterns   from    (110)    sputtering    (figures 
6l— 6i+)    tend    to    show    more   of    an   oval   outline    than   a   uniformly 
dense    central   oval   region.         The    central   region    can   be    made 
more   dense   by   including   atoms    sputtered   with   less    than    1    ev 
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perpendicular    energy   in   the   probability-of-sputtering   selection 
process.       This    did   not    seem   justifiable   for    (110)    pattern 
production    since    these    small    energy   atoms    were    excluded   from 
(100)    and    (111)    patterns.       The    1    keV    point   plot    (figure   6la) 
showed   indications    of    (100)    spots    but    these    cannot    be    clearly 
seen   in    the   deposit    pattern    (figure   6lb).       Their    presence    was 
substantiated  by   numerical   data.       The   average    energy   deposit 
data   showed    that    the   average    energy   of    atoms    in  ^100}  regions 
was   nearly   double    that   of    atoms   found   in    the    central    oval. 
This    is    a    constant    characteristic    of    non-central    spot   regions. 
At  higher    ion   energies    there   was   no    direct   indication    of    the 
presence    of    these    spots    either    by   numerical   data   or    energy 
deposit   data.      A   few    atoms   ha.ving  high    energy    were   found    in 
these   regions   but    the   area   density   was    too    small    to   form    a 
spot   in    the   pattern.       The    set   of    atoms    found    to   form    the 
central   oval   for    sputtering   at   all   ion   energies    was    not   a    well 
defined   group    such   as   has    been   found   for    (100)    and    (111) 
sputtering.       The   atoms   forming   the    (  100 )    spots   from    1    keV 
sputtering   are    those    which   are    either   n.n.'s    to    the    target 
or   n.n.'s    to    the   atom    which   initiates    sputtering   in   each 
horizontal   row. 

The    sputtering   ratio   for    (110)    sputtering   is    shown    in 
figure   65.      Agreement    with    experimental   data    from    th< 
numerical    and    curve    shape   aspect    was   poorer   for    this    surfac< 
than    the    other    two    surfaces.       This    is    apparent1  used    by 
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(110)    sensitivity   of    the    sputtering  ratio    to    the   A   parameter 
in    the    argon-copper    potential   function.       Levy  obtained 

better    curve    shapes    using  A    =    11.2+35    whereas   A   =    12.56 
has   been   used   in    this    study.       This    sensitivity   has    not  been 
observed   for    (100)    and    (ill)    sputtering. 

Energy   distributions    are    shown   in   figures    66-69.       The 
notable   feature    of    (110)    sputtering   is    that    the    energy   of 
sputtered   atoms    is    generally  higher   for    (110)    sputtering   than 
(100)    or    (111)    sputtering. 

D.     Results    of   Ancillary   Studies. 

During   the    course    of    the    simulation,    it   became   apparent 
that    certain   extensions    of    this    study    should  be    made.       It    was 
particularly   desirable    to   probe   five   areas: 

(1)  Sputtering   mechanisms    peculiar    to   high   bombard- 
ment   energies    (up    to   40    keV )  . 

(2)  Sputtering   of    the    (110)    surface  by   high   energy 
ions    which    are    directed    towards    (110)     channels. 

(3)  Sputtering    mechanisms    for    an   ion  heavier    than 
the    target. 

(4)  Sputtering     of    the    (0001)    basal    surface    of    a 
hexagonal    close   packed    crystal. 

(5)  Sputtering    the    (100)    surface    at    ev   ion    energies. 
Investigation    of    the   first    two    of    these   five   areas    were 

limited    somewhat   by    the    restrictions    for    the    copper— copper 
and   argon-copper   potential   function.       They   are   usually   assumed 
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valid   up    to    about    10    keV,    however,    if    small    impact    parameters 
are    assumed    to    rarely   occur,     valid    results    may    be    expected. 
Areas    (3)    and    (4)    involve    the   use   of   unknown   parameters. 
Accordingly,    the   validity   of    results   from    investigating   these 
two   areas    cannot   be   assured. 

The    search   for    mechanisms    peculiar    to   high   energy   sputtering 
was    made   using   the    (ill)    surface.       This    surface   was    considered 
the   one    most   likely   to    show    additional    mechanisms    since   it   is 
the    most   densely   packed   surface.       The   failure    to   find    mechanisms 
(part  A   of    this    section)    peculiar    to   high    energy   sputtering 
was   not    too    surprising. 

The    sputtering   of    the    (110)    surface   by   ions    directed   into 

(110)  channels    was    investigated  by   determining   the    sputtering 
ratio   for   ions   impacting   at   the    two   impact   points   located 
nearest,   the    channels.       The    sputtering   ratio    curve   for    channel 
shots    is    shown   in   figure    70.       The   discontinuity  between    the 
low    and   high   ion   energy    curves    is    assumed    to   be    the   result    of 
breakdown   of    the   potential  functions.       The   important   result- 
is    that   sputtering   always    occurred   for1    channel    shots,    even 
when  1+0   keV    ions    were   used. 

The    third   area    was    examined   using   xenon    to    sputter    the 

(111)  surface.      Argon-copper   potential   parameters    were   used 
since    those   for    xenon-copper    were   not   available.       Runs    were 
made   at  3    and    5   keV .      At    both   energies,    the   ion   pen<  i     ;ned 

at   least    to    the   fourth   layer    of    atoms;    no   indica tion   oi'    n.omti.rum 
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reversal    of    the   ion    was    observed.       The   ion   was   found    to   be 
channeled,    just   before    penetrating   the    third    layer,    into    (110) 
channels.       It   initiated    cascades    which    propagated   into    the 
crystal. 

The   fourth   area   was    investigated   to    see   if    the    effect    of 
surface    geometry    on    the    sputtering   deposit    patterns    could   be 
determined   by   a    comparison   of    (111)    sputtering   of   a   face 
centered    cubic    crystal   with    sputtering   the    (0001)    surface    of 
a   hexagonal    close-packed    crystal.      Hasiguti,    Hanada,    and 
YamaguchP J        have    sputtered    zinc    with   8   keV    argon;    the 
deposit    pattern    showed    an    outlined,    equilateral   hexagon    with  ■ 
a    central   haze.      An   attempt    was    made    to   reproduce    this 
pattern    with    the    simulation    model    modified   for    zine-zinc 
interactions.       The   potential   parameter   was    adjusted   for    a 
Born-Mayer    type    potential   V(r)    =  Ae"    '        where   A    =    52 
(23*23    )  keV.       This    relationship    was    determined   by  Ander- 

sen  and    Sigmund    -*  '    .       The    argon-zinc    potential    parameters 
(assuming   a    Born-Mayer    type    potential)    were    approximated   by 
using    those   for    argon-copper    interactions.       Neither    the 
sputtering   ratio    nor    deposit    pattern    matched    the    reported 
results.       This    was    not    unexpected    since    the    zinc    atoms    in 
the    crystal    must    be    represented   as    ellipsiods    (with    the    minor 
axis    in   the   basal   plane)    rather    than   represented   as    spheres. 

The    (  100 )    surface   was    sputtered   at    100    ev    to   investigate 
the    possibility    of    a    45    degree    pattern    rotation  for    low    ion 
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energy    sputtering.       The    number    of    atoms    which    were    sputtered 
was    extremely   small.      All    sputtered   atoms   had    energies    less 
than   10    ev.       In   was   found    that   the    atoms    sputtered    most 
frequently  from    the   regular    surface   were   n.n.n.'s    (figure    7i  )  • 
Bombarding   vacancy   and    stub    surface    conditions    also    resulted 
in   a   predominance    of   n.n.n.'s   being   sputtered   rather    than 
sputtering   of    the   atoms   forming   the    square   about    the    target 
atom . 

E.     Results    Common    to   All    Surfaces. 

Momentum   focusing    was    observed    to    occur    only   into    the 
crystal   for    the    three    copper    surfaces    studied.       This    was 
particularly    evident   from    atom    track   displays    o±     (111)    surface 
sputtering.       If    an   atom    underwent    momentum    reversal,    it   was 
always   by   reflection   from    atoms    located   no    deeper    in   the 
crystal    than    the    third   layer.       Even   for    '/    keV   ion   energies, 
crystal   atoms   located   deep   in   the    crystal   were   always    driven 
inward.      A    surprisingly   large    number    of    atoms    with    large 
energies   were   found    to    move    between   second   and    third   layers 
with    their    motion   nearly    parallel    to    the    surface    >;;nd    with   a 
small,    inward— directed   momentum   component  .      The   energy 
which   appeared   to   be   delivered   in  impulses    bo   atoms    still    a1 
lattice    sites,    dissipated    through    the    crystal.       Thi  s<       Ltoms 
rarely    acquired    more    than    kJU    ev    through    the    energy    impulse 
process.       Their    motion    was    restricted    since    the}-    were 
surrounded    by   other    atoms. 


The    target    atom    in    the    (100)    and    the    (110)    surfaces, 
when   struck   at   near    zero    impact   parameter   by    the   ion, 
always    transmitted    the    majority    of    its    energy   billiard-ball 
fashion   in   close   packed   rows    perpendicular    to    the    surface. 
(Atoms    to    which    this    energy    was    transferred    escaped    through 
the   back   face    of    the    crystal.)       There    were    a   few    cases    in 
which    one    of    these    atoms    would   initiate    another    chain.       The 
chain   was    never    more    than   2-3    atoms    in   length   before    the 
energy    was    either    wholly   dissipated    or    the    chain    stopped   by 
divergence    of    the    momentum    to   form    numerous    small    cascades 
directed   into    the    crystal   interior. 

The   unique   arrangement   of    atoms    in    the   first   layer   was 
observed    to   be    a   dominant   factor    in   determining   ejection 
directions   for   all    surfaces.      Atoms    which    were    ejected   other 
than   nearly   normal    to    the    surface   were    always    influenced   by 
their   neighbors.      Atoms    which    were    ejected   at   angles   near 
45    degrees    to    the    normal    would   reflect   from    their    neighbors, 
often    ending    up    as    normally    ejected   atoms.       This    effect    was 
most   pronounced   for    (111)    sputtering   since    the    six   atoms 
acted   as    a   lens,    but   it   was    also    seen   for    (110)    sputtering 
in   which    a    lens    is    formed   by    second    layer    atoms.       The    (100) 
surface    showed   a    strong   lens    effect   although    one    would   not 
necessarily   assume    that   it   occurs.       Sputtered   atoms, 
originally   adjacent   in    the    (100)    surface,    were   often   observed 
to   be    ejected   almost    simultaneously.       In    these    situations    they 
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were    ejected    nearly    normal    to    the    surface.       This    may    be    the 

(  OE.  ) 

genesis    of    dimers    recently    observed   by  Woodyard    JJ    . 

The   spot    regions    of    the   deposit   patterns   from    (100)    and 
(111)    sputtering   always    contained   more   high    energy    (^10    ev) 
than   low    energy   atoms.       But   the    energy   distribution    was    more 
uniform    than   increasing    with    the   ion    energy.       It    was   not 
peaked   at   any   one    or    group   of    energies.      Atoms    with   low    ion 
energies    were   predominantly   found   in    €he    central   region   of    the 
pattern. 
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6 .       Conclusions. 

The   observable    quantities    of    sputtering   appear    to   be 
interdependent     on   only   the    sputtering   mechanisms.       There   is 
no    evidence    that    these    quantities    are   interrelated.       The 
qualitative    and   quantitative   data   from    the    simulation   indicate 
that   the   deposit   pattern,    sputtering   ratio*       and   energy 
distribution   of    sputtered   atoms    cannot   be    correlated   with 
each   other;    one    cannot   predict   a    sputtering   ratio   from    an 
energy   distribution.       But,    each    of    these    quantities    can   be 
cross-correlated   between    surfaces. 

The    main  features    of    deposit   patterns    appear    to   be 
determined   only  by   the    surface   geometry   of    the    crystal.       The 
formation   of    (111)    pattern   spots    is    attributed   primarily   to 
assisted   focusing   by   the   hexagonal   lens.       The   predominance   of 
trigonal   rather    than   hexagonal    symmetry   is    considered   a 
natural   result   of    the   brief    (112)    channeling   observed   in   the 
mole      mechanism.       The    appearance    of    streaks   between    the 
spot   pairs    and    the    appearance    of    a   hexagonal    spot   pattern   at 
2   keV    suggests    that   use   of    a   hemispherical    collector    in   the 

simulation   will    show    the    three    {111+)    spots    which    were    observed 

(  qq  ) 
by   Southern   and   Robinsonv  ■J-J    .       However,    no    distinction    can 

be   made  between   the    (110)    and    (114)    spots    on   the   basis   of 

pattern   location.       If    the    mole  mechanism    is    predominant   at 

certain    energies    the    (110)    spots    should   be    more   intense    than 

the    (  11J+ )    spots.       If    the    squeeze    mechanism    is    predominant, 
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the    (110)    and    (114)    spots    should  have    more    equal   intensities. 
The    nearly    exact   simulation   of    a    (100)    pattern   at    3    keV 
is    considered   one    of    the   best   arguments    that    sputtering   is 
mostly   a    surface   phenomenon.       Neither    the   potential   form 
nor   its   parameters   for    copper-copper    and   argon-copper   inter- 
actions   are   known   with    certainty.       If    sputtering    was    a   deep 
phenomenon,    one    might   accept   an   argument    that    spots    in    the 
pattern    could  be   produced   without   exact   knowledge    of    the 
potential.       However,    the   hypocycloid    outline,    definitely   present 
in    the    simulation   pattern,    is    considered   the   feature    of    the 
pattern   which    would   be    most    sensitive    to    small   variations   in 
the   potential   if    sputtering   involved    more    than    the   first   few 
layers.       This    conclusion   is    consistent    with    the    apparent 
rotation   of    the    1    keV    pattern.      Atoms    found   in  ^2.1.1^  regions 
were    sputtered   when    the    energy    transferred   by   the   ion    to    the 

target   was   as    small   as    40    ev.       This    compares   favorably    with 

o 
a    calculated    transfer   value    of    75    ev   for    a    0.5   A     impact. 

parameter    \^^)  ,      A    sputtering    threshold    energy   of    50    ev 

has   been   reported    ^3«)# 

The   inability    to    produce    consistent    (100)    patterns    at   a 

given  binding   energy   for    5    and    7    keV    sputtering   is    undoubtedly 

due   to    the   unknown    condition   of    the    surface.       It    would   be 

unreasonable    to    assume    that    the   binding   energy   is   a    decreasing 

function   of    ion   energy   only   for    the    (100)    surface.       The 

disruption   of    the    surface    by    the   first    group    of    incident    ions 
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is    a   factor    which    cannot   be   ignored.       The    sputtering   ratios 
determined   for    the    (100)    surface   as    well   as   for    the    (110) 
and    (111)    surfaces,     were    within   reasonable    limits    when    the 
assumed   binding    energy   is    in    the    range    1.50    to    3-50    ev.       No 
definite    value    of    binding    energy    can   be    determined    from    the 
patterns    and   sputtering   ratios    unless    a   weighting  factor   is 
used   for    the    surface    condition    at   the    time    of    ion    impact. 
The    assignment    of    weighting   factors    would   be,    at   best,    a 
guess . 

Three    surprising   features    were    observed    in    simulation 
deposit   patterns.       The   first,    the   apparent   rotation   of    the 
(100)    pattern   at    1    keV ,    will   be   further    discussed.       The 
rotation   of    the    ( 100 )    pattern   for    low    energy   sputtering  has 
not   been    observed.       The    results    of    the    simulation    indicate 
that    such    a    rotation    is    possible.       This    conclusion   is   based    on 
empirical    rather    than    theoretical    considerations.       The   hypocycloid 
outline    found    at    3    keV    in    the    simulation    appears    to    be    formed 
by   the    same   atoms    which   form    (211)    spots    at    1    keV. 
Furthermore,    no    explanations   have   been   previously   proposed   as 
to    why   the    outline    is    a   hypocycloid    in    experimental    patterns. 
One    would    expect    that    a    circular    or    perhaps    a    non-distinct 
outline   would   be    observed   when    the    (100)    surface    is    sputtered. 
Accordingly,    it    is    proposed    that    the   hypocycloid    outline    is    a 
result    of    the    inability    of    the    (100)    surface    to    completely 

focus    atoms    sputtered   by   a    squeeze    mechanism    for   bombardment 
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energies    at    which    the    energy   propagation   is   not    confined    in    the 
surface. 

The    second   feature    was    the   unique   form    of    the    bonds    which 
connected    the    spots    in    the    simulation   pattern  for    (111) 
copper    sputtered   at    5    keV  .       The   appearance    of    the    similar 
outline,    a    spherical    triangle,    was    observed   in   a    pattern   of     (ill) 
copper    sputtered    with    1.5    keV    krypton    by    Yurasova    and 
Bukhanov' -^    .       The    third    feature    was    the   hexagon    spot   pattern 
at    2    keV,    a    pattern    seen    by   Anderson   and  Wehner*        '    for 
(111)    copper    sputtered    with    400    ev    mercury.       The    similarities 
observed   in    the    simulation   patterns    with    these   anomalous 
features    of    the    experimental   patterns    suggested    that    common 
factor    other    than    the    (111)    copper    surface    might    be   present. 
The    momentum    ratios    of    1+00    ev    mercury-2   keV         argon   and 
1. 5  keV    krypton-5   keV         argon   are    1.0    and   U.8    respectively. 
This    is    not    considered    to    definitely    establish   a    sputtering 
correlation   based   only   on   ion    momentum.       The    momentum    ratio; 
for    1.5    keV    krypton   ana   4    and   3   keV    argon   are   0.9    and    1.0 
respectively.      Ideally,    the    spherical    triangle    would   be   obst     \    d 
in    the   3    keV    (ill)    simulation   pattern.      Additionally!    wf-u 
(111)    copper   is    sputtered   at    high    temperatures,     the    triangu] 
outline   has    been   observed    to   become    more    pr    1    iiuiced    '     H      . 
The    absence    of    such    a    momentum    scaling    effed     1    >uld    bt 
easily    shown    by    experimentally    sput  tering    (111)    c  >pper    with 
ion    momentum    comparable    1    -    thai     oj     l+QQ    ev    mi  jr. 

6; 


The   general    consistency   of    the    numerical   values    of 
sputtering   ratios,    and    energy   distributions    of    sputtered   atoms 
and   likeness    of    the    patterns    with    experimental   data   is    a 
strong   argument   for    the    validity    of    a    computer    simulation    of 
sputtering.       It   is    remarkable    that    the    consistency   is    as 
good    as    it    is.       The    model    uses    only   a    repulsive    potential,    the 
crystal    size    used    in    the    simulation    is    an   infinitesimal    portion 
of    the    smallest   laboratory    specimens,    and    the    potential   form 
and   parameters    are    comparatively    crude.       It    is    concluded 
that    the   results    obtained    substantiate    the    concept   of 
transparency   and    the    occurrence   of    momentum   focusing   only 
within    the    surface    layer. 
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APPENDIX    A 
The    Beam    Model 
A.     Impact   areas    and   Impact   points. 

Each    crystal    surface    contains   an   intrinsic,    plane    geometric 
shape;    hexagon    for    (ill)    surface,    square   for    (100)    surface, 
and   rectangle    for    (110)    surface.       A    volume    element    of    the 
(hkl)    surface   is    defined  by   the   area   of    this    intrinsic    shape 
and   a   depth   of    some    number    of    (hkl)    planes.       This    volume    is 
chosen   so    that   by   translation   along   axes    of   a    Cartesian 
coordinate    system    the    entire    crystal    may  be   generated.       The 
smallest   intrinsic   area   which    can   be    chosen   for    each    surface 
and    still    satisfy   the    translation   requirement   is    shown   in 
figure    72.      A  finite   number    of    points    is    symmetrically   distributed 
within   each   of    these   areas    to   represent    the   infinite    set   of 
possible   points    of    impact  for   an   incident   ion.       These    smallest 
areas    are   further    divided   into    representative   impact   areas 
which    are    seen   from    figure    73    to    be    degenerate    under    appro- 
priate   rotations    and/or    inversions    of    the    coordinate   axes. 
Since    the    sectors    are    degenerate,    the    set    of    impact    points 
in   each    sector   is    also    degenerate;    only   the   points    contained 
in   one   impact   area   need   be    used    to    represent   bombardment   of 
the    entire    area.       The    coordinate    axes    rotation   and    inversion 
schemes    are    discussed   in  Appendix   C    (Deposit    Pattern 
Production)  . 

The   independence    of    the   impact   point    set   used,    with 
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respect   to    number    and   kinetic    energy   of    atoms    sputtered,    was 
tested   using    two    sets    of    points.       No    dependence    was    found. 
Results   for    the    (111)    surface   using   a   3   keV    ion  are    shown   in 
figure   74   as   an   example.       Set    1    points   are   those    shown  in 
figure    73 »    set   2    (not   shown)    is    a   set   of    eleven   points    which 
are   located   between    the   points    of    set    1. 
B.     Neutralization   of   Beam    ions. 

The   assumption    that:   argon   ions   are   neutralized   prior    to 

impact   on    copper    may    be   inferred   from    results    of    a    theoretical 

(   \1 ) 

study   of    secondary   electron   emission  by   Harrison    et.    al.  . 

Consideration   of    atom-atom    rather    than    Ion— atom    interactions 
gave   results    in   reasonable   agreement    with    experimental   data. 
A   supporting   argument,   for   neutralization   is   based    on 

/   -J  u  \ 

Hagstrurn's    theory   of   Auger    ejection   of    electrons  .       The 

probability   of   an   ion   being   neutralized   in   dx  at   x   is: 
Ft(x,v)»a   exp  J   -  exp     l-a(x-xm)        -a(x~xm)V» 
where    xm    ■    (l/a)  JLn     (A/av)    is    the    value    of    x    where    P+.    is 
maximum.       The   parameters   A   and   a   occur   in    the    transition 
rate  function,    and  _v   is    bhe    velocity   of   an  ion  starting  at   x      oo . 
Hagstrurn    used    tungsten   as   an   example   and   obtained    the    va)u< 
of   A   by   empirical    means,     the   value   of    a   from    published    data 
Rather   than   a1  bempl    an    exact   proof   of  neutralization    for   argon 
on   copper,    Hagstrurn's   results  for  argon   on    tungeti  n    have   been 
used   to   give  an  order   oi    magnitude,    at  worst,    approximation! 
Accordingly,    1.6j(  x1Tl^  r  .  I  tor    l     to    10    keV    argon    ion 
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copper;    this    is    a    reasonable    distance    from    the    crystal    surface 
and    one    may   assume    that    neutralization    occurs. 
C.      Ion    Deflection   by   Surface    Potential. 

As    a   first   approximation,    it   is    assumed    that   a   singly 
positive-charged   argon   ion   sees    the    crystal   as   a    single,    fixed 
scattering   center   having   2-1     positive    charge.       The    well   known 
equation   for    a    central   force-induced   hyperbolic    orbit   is    used 
in   conjunction   with   figure   75    to    determine    the    deflection. 

r(0)    =    (L2/mk)    /[-l      +    (H2EL2/mK2)l/2    COs(  Q->  &*    )\ 
L    =   ion   angular    momentum  m    =   ion    mass 

E    =    ion    energy  K   =    Q     Q    A^S© 

\J0  -   angle    of    closest   approach 
The    angular    momentum    L    is    determined    at    r-»co^^-*^" 
to   be   L    =    mvb   where    mv   /2    =    E    and   b    is    the    impact   parameter. 

The    equation   for    r(  0  )    is    put    in    a    more    convenient    form    by 

?  2  ? 

making-     the    substitution   L      =    (mvb)       =   2mEb      to    give 

r  (0)  =  b  (2Eb/k)/i  -1  t  [l  +  (2Eb/k)2l1/2  cos(<9~<90  ) 
The  angle  of  closest  approach  (Sfe  >  is  determined  as  a  function 
of    the    product    Eb     by    requiring    the    denominator    to    vanish    for 

(£>:r  ~ff     .       Once  Q0    ( as    a   function   of    Eb)    is   known,    the 
deflection    ratio    r/b      is    calculated   at    the    crystal    surface. 
The    deflection    calculated   is   for   an   unneutralized   ion,    thus 
greater    than    that   for    an    ion    which   is    neutralized   at    some 
distance   in   front    of    the    surface.       Plots    of    r/b    at    the    surface 
and  (0  as    functions    of    the    ion    kinetic    energy-impact    parameter 
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product   are    contained    in   figure    76.       The    numerical    values    are 

+  + 

for   an  Ar     —  Cu      system. 

A    more    sophisticated   approximation,    scattering  from    a 

fixed   dipole   of    same    charge    sign,    requires   a   messy   integration. 

Rather    than  follow    this    line,    one    may   expect   from    the   nature 

of    the   problem    that   the   path   of    the   approaching  ion   will   undergo 

some   oscillatory    motion   or   perhaps    corkscrew    motion   as    the   ion 

is   influenced   by   the    surface   potential.       In   either    case   one 

would   expect   that    the     net     acceleration   of    the   ion   parallel    to 

the    crystal   surface    would    be   no   greater    than    that   due    to   a 

single   fixed    scattering   center.       With    these    considerations, 

figure   76   is   used   to   determine   the  percentage   of   ions   which 

will  be  appreciably   deflected.      As   an   example   of   an  appreciable 

ion  deflection,    consider   r/b   ■   2.0  |    from   figure   7b,    Eb   «   3  •.  J5    * 

o 
10"'    keV-A.      An   energy   range   of    1-10    keV    for   argon   ions, 

that   used  in  this   study,    corresponds   to   a   range   of   Impact 

o 
parameters   of   Q.Q35   -   0.0035  A.      'The   fraction   oi    surface 

area,    and   therefore   fraction   of    ions    which   will   be   appreciably 

deflected  is    (0.035)    /( 1.26)      -  7.7   x   10~Jf  for   1   keV    Ions   .ma 

o 
7  >7   x  10~°  for   10   keV    ion©    (1.2b  A  is   Che   copper   atom 

effective  radius   In   the   crystal).      This   negligibly   small   fraction 

of  ions   cannot:   influence    the   macro  aspects   of    the   spul  terij 

model. 

D.     Equilibrium   fc3tate   During   Impact. 

The  use   of   a  single  atom   approximation    to    the   beam    wan 
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the    crystal    model   described   in    the    main    text   requires    that    the 
crystal   region  be    in   an   equilibrium    state   at   the    time    of    an 
ion   impact.       It    is   not   necessary   that   this    equilibrium    state 
be   identical    to   previous    equilibrium    states    since    the   variety 
of    surface    configurations    available    in   the    model   provides   for 
random    surface    conditions.       It    is    necessary   that    the    time 
required   for    the    crystal    region    to    return    to    an    equilibrium 
state   be    small   with    respect    to    the   arrival-time   intervals    of 
the   ions.       Satisfaction   of    this    condition   is    determined  by 
comparing   an   experimental  beam    flux   over    the   area   of    the 
crystal   face    of    the    model   used    to    the    time   required   for 
completion   of    all    energetic    collisions    in   the    model. 

Beam    intensities    of    the    order    of    lOOwamps/cm      were   used 

(IP) 
by   Magnuson   and    Carlstonv  x    '  .       If    a   beam    of    this    intensity 

o 

is   incident   on  a    crystal    surface   area   of    less    than    lO-^A^    such 

as   in   the    model,    the   ion   flux   over    this    area   is    less    than    100 

ions/sec,    an   ion   arrival— time   interval   of    10—      seconds;    all 

energetic    collisions    in    the    model   are    completed   within  about 

—12 
10  seconds.       The   relaxation    time    of    the    region   is    thus    much 

smaller    than    the    ion   arrival-time   intervals    and    the    condition, 

that   the    crystal   region   be    in   an   equilibrium    state   at   the    time 

of   impact,    is    satisfied. 
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APPENDIX    B 

Positioning    the    Ion 

The    ion    is    positioned    tangent    to    the   first    target    atom    it 

will    strike.       Since    this    is    a   dynamic   rather    than    static    process 

it   is    not   necessary   that    this   be   a    stable   position   on    the    crystal 

surface.       Figure    77    shows    the   ion   at   its    arbitrary   initial 

position   and   calculated   final   position,    both   -with   respect    to    the 

impact   point   and    target   atom.       The    initial   position   is   a    small 

distance   in  front   of    the    surface,    beyond    the   eroded   potential 

range    of    the    crystal   atoms.       This   position   is    described   by   a 

vector    r-     originating   at    the    impact   point   and  having   direction 

parallel   but   opposite    to    the   ions    velocity   vector.       The    vector 

rg>    from    the   impact   point    to    the    target   atom,    is   known 

since    the    target   atom's    coordinates    are   known;    the    vector   r1 

is   known,    and    it   is    desired    that   vector   r        have   a    magnitude 

equal    to    the    distance   between    centers    of    two    tangent   atoms, 

2ra.     The    law    of    cosines    is    used    to    compute    the    magnitude    of 

r~    which    lies    along   r    .      Accordingly,    the   following   calculations 

are   made: 

cos  Ot  -       r       .    r 
1  2 


I'll   M 

(r-a)      -    (2r_    cos*    )      (r,  )     \    (r      -    (2r    )*")    =   0 
J  «  3  2  o 

(xvJ    ••    c2    cosoi     i         r|    cos2*.    -    (r2    -    (2r    )    ) 
The  positive    square    rool     solution  is    chosen    Lu   give    the    tangent 
on    ili;   outsid<     hemisphere   of    the    target    atom. 
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The   ion   positioning   is    accomplished   automatically   for    each 
run  by   subroutine    START,    which   is    contained   in   the    computer 
program    for    the    sputtering    simulation. 
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APPENDIX    C 

Production   and  Analysis    of    Sputtering   Deposir    Patterns 
A.      Production. 

A    sputtering    deposit:   pattern   represents    the    intersection 
points    of    atoms'    velocity   vectors    with    the    surface    of    a 
collector    plate.       The    sputtered    atoms    from    the    simulation 
are    collected     on    a   flat    plate   by   determining    these    points    of 
intersection.       Each    atom    -which    exits    through    the    crystal 
surface    is    initially    assumed    to    have   been    sputtered    and    a    data 
card  has    been   prepared   for    each   one.       Data    of    particular 
interest   are    each    velocity    component    magnitude   and    the   kinetic 
energy   perpendicular    to    the    surface;    data    of    secondary 
interest   are'  atom    number,    impact   point   used   and   ion   kinetic 
energy.       (The    use    of    these    last    data   will  be   discussed    shortly.) 

It   is    recognized    that    each    sputtered   atom   has    lost    some 
energy    to    overcome    the    surface   binding   energy.      An   assumed 
value    of    binding    energy   is    subtracted   from    the    perpendicular 
kinetic    energy   and   a   new    perpendicular    velocity    component   is 
calculated.       Parallel   velocity    components    are    normalized    to    the 
new    perpendicular    component    to   give    a    two    dimensional    coordi 
ate   point.       This    point    is    the    intersection    of    the    atom's 
velocity   vector    with    an   imaginary    collector   plate   placed    at    ur;ii 
distance   from    the    tar-get    surface.       Each   point   is    then    ro1        >-'d 
and/or    mirrored   about    the    coordinate    axes    to    give    the    inter- 
section  of    a    velocity    vector    which    would  have    resulted    for    an 
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ion   impacting  at    the    corresponding   impact   point   in   each    of 
the    other    impact   areas.       The   impact   areas    are    shown   in 
figure    73    and    the    coordinate   point   rotation   and    mirroring  values 
are    listed    in   program    DATASORT    which    is    used    to    generate 
the    points.       Each    point    is    plotted    using   program    DATAPLOT 
with    a    CDC160A    computer    and    a    CalComp    plotter    to    give    a 
point   plot.       The    dimensions    of    the    plots    are    in    target-to- 
collector    or    T-Q    units    since    these    points    have   been   normalized 
to    unit    target-to-collector    distance.       The    scale   which   has 
been   used   permits    plotting   deposit   points    of    atoms    which   have 
been   sputtered   within   an    escape    cone    of    about   63    degrees 
(57      degree    cone    shown    in   figures).       This    has    been   found 
satisfactory   to    contain   all   pattern  features    of    interest. 

The   point   plots    show    only   point   patterns    and   therefore 
do    not   accurately    simulate    experimental   patterns    which   are 
area   density   patterns.       The    conversion   of    a   point   plot   to    a 
smooth    area    density    pattern    is    made    by   photographing    the    point 
plot    with    the    camera   defocused    such    that   no    single   point    is 
distinguishable   but   high   and   low    density    areas    are   prominent. 
Developing   and   printing   is    controlled    to   bring   out    the  high 
density   areas    while    maintaining   the  haze   background.       Loss    of 
intensity    in    some    spot    regions    cannot    be    avoided    such    as    seen 
in   figures    6b    and    12b.       The    process    must   be    adjusted   for 
each    pattern;    the    sequence    of    photography,    developing,    and 
the    printing   of    the    positive    image    is    highly    dependent    on    the 
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ability   of    the   photographer.       Concurrently,     the    quality   of 
reproduction   of    these   patterns    in   printed   form    depends    on    th< 
plate   preparation,    the   printing  press,    and    the   paper    used. 
B.     Analysis. 

A   30    x  30    square   grid   is   placed    over    the    central    3.0    x    3*0 
unit   square   of    the   raw    pattern   by   program    DATAGRID. 
The   identification   number    of    an   atom    and    the   impact   point-    and 
impact   area   of    the   run   in   which    that   atom    was    sputtered   are 
recorded  for    each   atom    in    the   grid   square.      The    total   number 
of   atoms,    the    total   energy,    and   the   average    energy   per   atom 
for   each   grid   square   are   printed   in    separate   square   arrays. 
The   individual   grid    square   data   provides   for    correlation 
between   an  atom's    crystal   location   and   its    deposit   point   in 
the   pattern;    the   square   arrays    of    the   number    of    atoms    and 
the   total   and   average   energy   densities    may   be    compared 
directly   with    the   point   plot   or    smooth   pattern   for   analysis    of 
pattern  features. 
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APPENDIX    D 
2.       Section    1:       Memory    block    allocations    and   functions    are 

established.       All    storage    cells    are    zeroed. 

Input   and   output   formats    are    specified. 

Constants    are    set,    and    the    target    material, 
ion    species,    and    crystal   face    to    be    sputtered 
are    read    in. 
Section   2:       The    ion    energy,    impact    area,    and    impact 

point    are    read    in.       Constants    peculiar    to    the 
run    are    established.       The    appropriate    sub- 
routine  is    called    to    build    the    crystal,    and 
crystal    boundaries    are    calculated.       Subroutine 

START    is    called    to    position    the    ion.       The 
initial    value    of  At    is    calculated.       Initial 
coordinates    of    all    atoms    are    assigned,    and 
their    velocity    components    are    zeroed    (  except 

for    the    ion).       The    initial    coordinates    of    the 
crystal    atoms    are    printed. 
Section   3*       Forces    are    calculated   by    calling    subroutine 

STEP.       Atoms    are    moved    to    their    intermediate 
positions    in    the    two    step    cycle.       Subroutine 

STEP    is    called    again,    and    atoms    are    moved 
to    their   final   positions.       Final    velocities    i  re 
then    computed.       Force    components    are 
zeroed   m   preparation   for    the    next     timestep. 
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The    maximum    kinetic    energy   is    determined  for 
the    calculation   of  b*  t . 
Section   4*       The    time    remaining   until    cutoff    is    determined. 

If    there   is    insufficient    time    to    complete 

another    timestep,    terminal   data   is   printed. 

Potential    energy   is    calculated   and   summed 

with   kinetic    energy    to   give    the    total    energy 
for    energy   balance    check    (manual    check).       The 
data   for    atoms   having   potential    energy  greater 
than    the    thermal    energy   are   printed. 
Section    5:      Atoms    which   have   kinetic    energy   but    not 

potential    energy   are   assumed    to   be   free   of 

the    crystal.       They   are   assigned   LCUT=1. 
The    surface   through    which   an   atom    exited   is 

determined   and   a    code   assigned.      A    maximum 
potential    energy   is   found    (for  Lt    calculation) 
among   atoms  which    do    not   have   LCUT=1.       If 
the    maximum    potential    energy   is    less    than   a 

minimum    value    the    terminating   process   begins. 

If    this    energy   is   greater    than    the    minimum 
value,    another    timestep   begins. 
Section   6:       Pertinent   data   for    all   atoms    is    printed.      Atoms 

which   have    exited    through    the   front    of    the 
crystal    or    will    exit    through    the   front    are 

assumed    to   be    sputtered.       Data   for    these 
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atoms  areprunted  and  a  data  card  punched 
for  each  atom.   A  new  data  card  is  then 
read  into  initiate  another  run. 
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APPENDIX  E 


AC 

AIX 
AIY   . 

AI2 

ALFA 

ALFA2 

B 

BENGY 

BULLET 

BX 

BZ 

CELS 

COX 
COY 
COS 

COXI 
COYI 

COZI 

CVB  ' 

CVD 
CVE 
CVM 

CVR 

DCOX 
DCOY 

DCOZ 

DFF 

DIST 


Glossary   for    FCCSPUT 
Parameter   for    target   force   function    correction. 

Floating   point   form    of    IX,    IY,    12. 

Cosine    of    the   angle  between   vectors    Rl,    R2. 

ALFA    squared. 

Reciprocal   of    magnitude    of    atom   velocity. 

Energy    which   an   atom    within    the    crystal   at    shut- 
down   must   have    to   be    considered    sputtered. 

Variable   representing   primary    material. 

Unsealed   x,    z    coordinates    of    the    impact   point 

Frictional   force    multiplier.       (See    CVS) 

Direction   cosines    of    primary   velocity   vector. 

Negative    values    of    COX,    COY,    COZ. 

A    constant. 

Converts    meters    to    angstrom    units. 
Converts    electron   volts    to   joules. 
Converts    atomic    mass   units    to    kilograms. 

Converts    lattice   units    to   angstrom    units. 

Direction    cosines    of    sputtered   atom    velocity   vector. 


Distance    difference   between   nearest    neighbor 
distance   and   actual   atom  ^differences. 
Distance   between    any    two    atoms. 
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DRX 
DRY 

DRZ 

DT 
DTI 


DTOD 

DTOM 

DTOMB 


DX 

DY 

D2 

ECUT 


EI 

EMAX 

EV 

EVR 

EXA 
EXB 

FA 

FAC 

FM 

FOD 

FORCE 


x,    y,     z    components    of    DIST. 


Length    of    timestep    in    seconds. 

Number    of    lattice    units    most    energetic    atom 
may    move    in    one    timestep. 

DT/CVD    -    a    ratio    used    to    avoid    repeated    division 

DT/PTMAS    -    a    ratio    used    to    avoid    repeated 
division. 

DT/PGMAS    -    a    ratio    used    to    avoid    repeated 
division . 


x,    y,     z    distances    atom    has    moved   from    initial 
position. 

A    lower    limit    on    an   atom's    potential    energy.       If 
energy    is    less    than    or    equal    to    ECUT    the    program 
shuts    down. 

A    cutoff    energy. 

The    maximum    energy    encountered    in    any    cycle. 

Primary    energy   in    electron    volts. 

Primary    energy    in   kilo— electron   volts. 

Potential   function    parameters. 

The    component    force    increment    on    an    atom. 

The    minimum    distance    the    primary    is    positioned 

in   front    of    the   first    xz    plane    at    start    time. 

A    small    number    used    in    checking   potential    energy 
zero    point. 

FORCE/DIST    -   a   ratio   used    to    avoid   repeated 

division. 

Numerical    value    of    the    force   function    with    a 
variable   parameter. 


1  I )  M 


FORF 
FPTC 
FFTF 
FRO 

FX 
FY 

F2 

FXA 

GMAS 
HDTOD 

HDTOM 

HDTOMB 

HGMAS 
HTMAS 
ICUT 


IHB 
IHS 
IHT 
IH1 

IN 


INDEX 


I  NOW 


Target   atom   force   function. 

The    corrective   force   value   at   ROE. 

The    corrective   force   function. 

Numerical   value   of    the   target   force   function   at 
ROE. 

x,    y,    z    components    of    total   force   on   an   atom. 

Force   function   parameter. 

Target   atom    mass    (in   a.m.u.) 

1/2    DTOD    -   a   ratio   used    to    avoid   repeated 
division. 

1/2    DTOM   -   a   ratio   used    to   avoid   repeated 
division. 

1/2    DTOMB    -   a   ratio    used    to   avoid   repeated 
division. 

1/2    GMAS    -   a   ratio    used    to    avoid   repeated 
division. 

1/2    TMAS    -   a   ratio   used    to   avoid   repeated 
division. 

Used    to    provide    output    prior    to    time    limit 
shutdown. 


IDUM  Dummy   variable. 


Alpha-numeric    arrays    for    titling. 


Odd— even   integer    used    to    determine   atom    site 
establishment . 

Integer    (0    or    1)    used   in   determining   dynamics 
cycle    step. 

Time   program   has   been   running   in    seconds. 
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IP  Subscript    value    of    atom.       Used      in        subroutine 

STEP. 

ISHUT  Time    left    prior    to    time    limit. 

ISPX  Fixed    point    values    of    SPX,    SPZ, 

ISP2 

ISS  Subscript    value    of    most    energetic    atom. 

IT  Unsealed   fixed    point    x    coordinate    used    in   lattice 

generation.       Also    a   dummy   variable    in   function 

TIMEF. 

ITT  Odd-even   integer    used    to    determine    atom    site 

establishment. 

IX 

IY  Number    of    x,    y,     z    planes    of    crystal. 

IS 

IXP 

IYP  Crystal   dimensions    in   x,    y,    z. 

IZP 

JSHUT  Cutoff    variable   based    on    total    potential    energy 

of    crystal. 

JT  Unsealed   y   coordinate   used   in    crystal   generation. 

UTS  Variables    used    to    establish   atom    sites. 

UTT 

KCUT  Identifies    exit    point    of    atom. 

KT  Unsealed    z    coordinate    used    to    establish    atom    site 

LCUT  Used    to    identify    atoms    which    are    not    included    in 

calculations . 

LL  The   highest   numbered   atom    in    the    crystal. 

LS  Sum    of    the    Miller    index    integers. 

LSS  Used    to    identify    type    of    surface,    i.e.,    regular, 

stub,    vacancy. 

M  An    integer    used    to    begin    atom    numbering. 


I  1 1 ) 


N  Subscript   of    the   atom    to   be   removed   for 
vacancy    surface. 

ND  Data    output   increment. 

NPAGE  Page    numbering   variable. 

NS  Initial    print    statement    cycle  . 

NSHUT  Cutoff    variable   based   on    too    long   a    timestep. 

NT  Timestep. 

NTT  Timestep   limit   before    shutdown. 

PAC  Same   as   AC    except   applicable    to    primary. 

PEXA  Primary   force   function   parameters. 
PEXB 

PFIV  A    constant   =   0.5. 

PFORF  Primary   force   function. 

PFPTC  PFPTF    evaluated   at    ROE. 

PFPTF  Primary    corrective   force   function. 

PFRC  PFORF    evaluated    at    ROE. 

PFXA  Primary   force   function   parameter. 

PGMAS  Primary    mass    in    kilograms. 

PKE  Kinetic    energy    of    an   atom. 

PKEY  Y    component    of    kinetic    energy    of    an   atom. 

PLA  Crystal    plane    (alphanumeric    variable). 

PLANE  Same    as    PLA. 

PNUM  Impact    point    (alphanumeric    variable). 

POT  Potential   energy   between   two   atoms. 

POTF  Target    potential   function. 

PPE  Potential    energy    of    an    atom. 
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RBX 

RBZ 

ROE 

ROEM 

PPOTF  Primary    potential    function. 

PPTC  PPOTF    evaluated    at    ROE. 

PRI  Chemical    symbol   for    primary    material. 

PTC  POTF    evaluated    at    ROE. 

PTE  Total    energy    of    an    atom    (potential    t    kinetic). 

PTMAS  Target    mass    in    kilograms. 

QM  A    small    number    used    in    checking   kinetic    energy 

zero    point. 

QUIT  Cutoff    variable    checked    against    total    potential 

energy. 

Unsealed    x,     z    coordinates    of    impact   area 
reference    point 

Nearest    neighbor    distance. 

ROE    -    DTI    (one    timestep    distance    less    than 
n.  n.    distance  )  . 

ROE2  ROE    squared. 

RX 

RY  x,    y,     z    coordinates    of    atom    at    any    time. 

RZ 

RXBND 

RYBND  x,    y,     z    coordinates    of    crystal    boundaries    other 

RZBND  than    zero. 

RXI 

RYI  x,    y,     z    coordinates    of    atom's    initial    position. 

R2I 

RXK 

RYK  x,    y,    z    coordinates    of    temporary   position    of 

R2K  atom    during    force    cycle. 

RXS 

RYS  x,    y,     z    coordinates    of    impact    point. 

RZS 

Rl  Vector    from    impact   point    to    Initial    primary   posit  i 
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R1X 
R1Y 
R1S 

R1R2 

R2 

R2SQ 
R3 

SAVE 

SCX 
SCY 

scz 

SLOW 

SPX 

SP2 

SSCZ 

SUR 
TAR 

TARGET 
TE 

TEMP 

THERM 

TI 

TIME 

TIMO 

TMAS 


x,    y,    z    coordinates    of    initial   primary   position. 

Scalar   product   of   vectors    Rl,    R2. 

Magnitude    of    vector    from    impact   point    to    first 
atom    hit   by    primary. 

R2    squared. 

Magnitude    of    vector    from    impact    point    of 
primary    start    position. 

1/2    POT. 

x,    y,     z    coordinate    scale    factors. 

Cutoff    variable    checked   against   a   long    DT . 

x,    z    distance   from    impact   area   reference   point 
to   impact   point. 

A    z    scale   factor   used   for  (ill)  plane    lattice 
generation. 

Plane    (alphanumeric   variable). 

Chemical    symbol   for    target    material. 

Target    material    (alphanumeric   variable). 

Total    energy    of    crystal    atoms    (kinetic    + 
potential)  . 

Temperature    of    lattice    in   degrees    Kelvin. 

Thermal    energy    of    atom. 

Computer    time    program    has    been    running. 

Elapsed    problem    time. 

A   function   to    convert    seconds    to    minutes. 

Target    atom    mass    in   kilograms. 
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TPKE  Total    kinetic    energy    of    crystal    atoms. 

TPOT  Total    potential    energy    of    crystal    atoms. 

VOL  Magnitude    of    primary   velocity   vector. 

VSS  Storage    variable    for    velocity    components. 

VX 

VY  x,    y,     z    components    of    atom    velocity. 

VZ 

X  Unsealed    x    coordinate    used    in    crystal    generation. 

XI T  Floating   point   form    of    IT    used    in   function    TIMO 

Y  Unsealed    y    coordinate    used    in    crystal    generation. 

2  Unsealed    z    coordinate    used    in    crystal    generation. 

ZE  A    constant    =    0.0 

ZP  Floating   point   form    of    JTT. 
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Erosion  of  the  force  function 

at  r  ^  2rQ    results  in  a  loss 

of  energy  in  the  model  .     This 

is  corrected  in  the  STEP  subroutine, 


Figure  1 
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Argon   Copper     Sputtering 
(  III )      Surface 
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Figure  6a 
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Figure  6b 
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Figure  7a 
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